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Abstract 

Traditional ideas for testing unification involve searching for the decay of the proton and its branching 
modes. We point out that several astrophysical experiments arc now reaching sensitivities that allow them 
to explore supersymmetric unified theories. In these theories the electroweak-mass dark matter particle 
can decay, just like the proton, through dimension six operators with lifetime ~ 10^^ sec. Interestingly, this 
timescale is now being investigated in several experiments including ATIC, PAMELA, HESS, and Fermi. 
Positive evidence for such decays may be opening our first direct window to physics at the supersymmetric 
unification scale of Mqut ~ 10^^ GeV, as well as the TeV scale. Moreover, in the same supersymmetric 
unified theories, dimension five operators can lead a weak-scale superparticle to decay with a lifetime of 
~ 100 sec. Such decays are recorded by a change in the primordial light element abundances and may well 
explain the present discord between the measured Li abundances and standard big bang nucleosynthesis, 
opening another window to unification. These theories make concrete predictions for the spectrum and 
signatures at the LHC as well as Fermi. 
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I. LONG LIFETIMES FROM THE UNIFICATION SCALE AND ASTROPHYSICAL 
SIGNALS 

One of the most interesting lessons of our times is the evidence for a new fundamental scale 
in nature, the grand unification (GUT) scale near Mqut ~ 10^^ GeV, at which, in the presence of 
superparticles, the gauge forces unify [1, 2]. This is a precise, quantitative (few percent) concor- 
dance between theory and experiment and one of the compelling indications for physics beyond the 
Standard Model. Together with neutrino masses [3, 4], it provides independent evidence for new 
physics near 10^^ GeV, significantly below the Planck mass of Mpi w 10^^ GeV. The LHC may con- 
siderably strengthen the evidence for grand unification if it discovers superparticles. Furthermore, 
future proton decay experiments may provide direct evidence for physics at Mqut- In this paper 
we consider frameworks in which GUT scale physics is probed by cosmological and astrophysical 
observations. 

In grand unified theories the proton can decay because the global baryon-number symmetry of 
the low energy Standard Model is broken by GUT scale physics. Indeed, only local symmetries can 
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guarantee that a particle remains exactly stable since global symmetries are generically broken in 
fundamental theories. Just as the proton is long-lived but may ultimately decay, other particles, 
for example the dark matter, may decay with long lifetimes. If a TeV mass dark matter particle 
decays via GUT suppressed dimension 6 operators, its lifetime would be 



Similarly a long-lived particle decaying through dimension 5 GUT suppressed operators has a life- 
time 



Both of these timcscalcs have potentially observable consequences. The dimension 6 decays cause 
a small fraction of the dark matter to decay today, producing potentially observable high energy 
cosmic rays. The dimension 5 decays happen during Big Bang Nucleosynthesis (BBN) and can leave 
their imprint on the light element abundances. There is, of course, uncertainty in these predictions 
for the lifetimes because the physics at the GUT scale is not known. 

If the dark matter decays through dimension 6 GUT suppressed operators with a hfetime as 
in Eqn. 1 it can produce high energy photons, electrons and positrons, antiprotons, or neutrinos. 
Interestingly, the lifetime of order 10^^ s leads to fluxes in the range that is being explored by a 
variety of current experiments such as HESS, MAGIC, VERITAS, WHIPPLE, EGRET, WMAP, 
HEAT, PAMELA, ATIC, PPB-BETS, SuperK, AMANDA, Prejus, and upcoming experiments such 
as the Fermi (GLAST) gamma ray space telescope the Planck sateUite, and IceCube, as shown in 
Table I. This is an intriguing coincidence, presented in section 2, that may allow these experiments 
to probe physics at the GUT scale, much as the decay of the proton and a study of its branching 
ratios would. Possible hints for excesses in some of these experiments may have already started us 
on such an exciting path. 

GUT scale physics can also manifest itself in astrophysical observations by leaving its imprint 
on the abundances of hght elements created during BBN. For example neutrons from the decay of 
a heavy particle create hot tracks in the surrounding plasma in which additional nucleosynthesis 
occurs. In particular, these energetic neutrons impinge on nuclei and energize them, causing a 
cascade of reactions. This most strongly affects the abundances of the rare elements produced 
during BBN, especially ^Li and '^Li and possibly ^Be. 




(1) 
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TABLE L The lower limit on the lifetime of a dark matter particle with mass in the range 10 GeV < < 
10 TeV, decaying to the products listed in the left column. The experiment and the observed particle being 
used to set the limit are listed in the top row. HESS limits only apply for > 400 GeV and are shown 
for two choices of halo profiles: the Kravtsov (K) and the NFW. PAMELA limits are most accurate in the 
range 100 GeV < < 1 TeV. All the limits are only approximate. Generally conservative assumptions 
were made and there are many details and caveats as described in Section IL 

In fact, measurements of both isotopes of Li suggest a discrepancy from the predictions of 
standard BBN (sBBN). The observed ^Li abundance of ^ ~ (1 - 2) x 10-^° is a factor of several 
below the sBBN prediction of ^ f=:i (5.2 ± 0.7) x 10"^° [5]. In contrast, observations indicate a 
primordial ^Li abundance over an order of magnitude above the sBBN prediction. The Lithium 
abundances are measured in a sample of low-metallicity stars. The Li isotopic ratio in all these stars 
is similar to that in the lowest metallicity star in the sample: ^ = 0.046 ±0.022 [6]. This implies a 
primordial ^Li abundance in the range ^ ?a (2 — 10) x 10"^^, while sBBN predicts ^ ^ 10^^'^ [7, 8]. 
The apparent presence of a Spite plateau in the abundances of both ^Li and ^Li as a function of stellar 
temperature and metallicity is an indication that the measured abundances are indeed primordial. 
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Of course, though there is no proven astrophysical solution, either or both of these anomahes could 
be due to astrophysics and not new particle physics. Nevertheless, the Li problems are suggestive 
of new physics because a new source of energy deposition during BBN naturally tends to destroy 
^Li and produce ^Li, a nontrivial qualitative condition that many single astrophysical solutions do 
not satisfy. Further, energy deposition, for example due either to decays or annihilations, during 
BBN most significantly affects the Li abundances, not the other light element abundances, making 
the Li isotopes the most sensitive probes of new physics during BBN. Finally, a long-lived particle 
decaying with a ~ 1000 s lifetime can naturally destroy the correct amount of ^Li and produce the 
correct amount of ^Li without significantly altering the abundances of the other light elements ^H, 
=^He, and ^Hc. 

In this paper we explore astrophysical signals of GUT scale physics in the framework of super- 
symmetric unified theories (often referred to as SUSY GUTs or supersymmetric standard models 
(SSM) [2]), as manifested by particle decays via dimension 5 and 6 GUT suppressed operators. 
In order to preserve the success of gauge coupling unification, we work with the minimal particle 
content of the MSSM with additional singlets or complete SU(5) multiplets. The effects of GUT 
physics on a low energy experiment are generally suppressed. However, these effects can accumulate 
and lead to vastly different physics over long times depending on the details of the higher dimension 
operators generated by GUT physics. For example, particles that would have been stable in the 
absence of GUT scale physics can decay with very different lifetimes and decay modes depending 
on the particular GUT physics. Conversely, such decays are sensitive diagnostics of the physics 
at the GUT scale. The details of these decays depend both upon the physics at the GUT scale 
and the low energy MSSM spectrum of the theory. So, astrophysical observations of such decays, 
in conjunction with independent measurements of the low energy MSSM spectrum at the LHC, 
would open a window to the GUT scale - just as proton decay would. The role of the water and 
photomultipliers that register a proton decay event are now replaced by the universe and either 
the modified Li abundance, or the excess cosmic rays that may be detected in todays plethora of 
experiments. 

To characterize the varieties of GUT physics that can give rise to decays of would-be-stable 

particles we enumerate the possible dimension 5 and 6 operators, each one of which defines a 
separate general class of theories that breaks a selection rule or conservation law that would have 
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stabilized the particle in question. This approach has the advantage of being far more general 
than a concrete theory and encompasses all that can be known from the limited low-energy physics 
experiments available to us. In other words, all measurable consequences depend on the form of the 

operator and not on the detailed microphysics at the unification scale (e.g. the GUT mass particles) 
that give rise to it. 

The presence of supersymmetry, together with gauge symmetries and Poincare invariance, and 
the simplicity of the near-MSSM particle content greatly reduces the number of possible higher 
dimension operators of dimension 5 and 6 and allows for the methodic enumeration of the operators 
and the decays they cause. This is what we do in Sections III, IV, and V. We work in an SU(5) 
framework because any of the SU(5) invariant operators wc consider can be embedded into invariant 
operators of any larger GUT gauge group so long as it contains SU(5). Of course, such operators 
may in general contain several SU(5) operators so the detailed conclusions can be affected. As an 
example, we study a model that is made particularly simple in an SO(IO) GUT in Section V. 

The dimension 6 operators of Section III have many potentially observable astrophysical signals 
at experiments searching for gamma rays, positrons, antiprotons or neutrinos. We separate them 
into R-conserving and R-breaking classes and in each case we also build UV models which give rise 
to these operators. The dimension 5 operators of Section IV have potentially observable effects on 
BBN, may solve the Lithium problems, and give dramatic out-of-time decays in the LHC detectors. 
Section V presents frameworks and simple theories that have both dimension 5 and 6 operators 
destabihzing particles to lifetimes of both 100 sec and 10^^ sec to explain Li and lead to astrophysical 
signatures today in PAMELA/ATIC and other observatories. 

In Section VI we similarly look at the consequences for Fermi/Glast and PAMELA/ATIC. In 
Section VII wc outline LHC-observable consequences of some of these scenarios (operators) that 
could lead to their laboratory confirmation. 

Finally, there are two more classes of theories that fit into the elegant framework of supersym- 
metric unification. The first of these are SUSY theories that solve the strong CP problem with an 
axion. The other is the Split SUSY framework. Both these frameworks provide particles that can 
have lifetimes long enough to account for the primordial lithium discrepancies without additional 
inputs, and are included in Section IV. 
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II. ASTROPHYSICAL LIMITS ON DECAYING DARK MATTER 



In this section we describe the existing astrophysical hmits on decaying dark matter, as 
summarized in Table I. Note that the hmits on dark matter decaying into many different fi- 
nal states (e.g. photons, leptons, quarks, or neutrinos) are similar even though they arise from 
different experiments. These different observations are all sensitive to lifetimes in the range 
given by a dimension 6 decay operator, as in Eqn. (1). This can be understood, at least for 
the satellite and balloon experiments, because these all generally have similar acceptances of 
~ (1 m^)(l yr)(l sr) 3 x 10^^ cm^ s sr. For comparison, the number of incident particles 
from decaying dark matter is ~ ^^'^ ^(0.3 ^^^^3 )(10~^^ s~^) ^ 10"-^° cm~^ s~^ sr~-^, where 
these could be photons, positrons or antiprotons for example, depending on what is produced in the 
decay. This implies such experiments observe ~ (3 x 10^^ cm^ s sr) x (10"-'^° cm~^ s~^ sr"-*^) fa 30 
events which is in the right range to observe dark matter decaying with a dimension 6 decay lifetime. 

Not only are there limits from astrophysical observations, there may be indications of dark 
matter decaying or annihilating from recent experiments. PAMELA [9] has observed a rise in the 
positron fraction of cosmic rays around 50 GeV. ATIC [10] and PPB-BETS [11] have observed a 
bump in the spectrum of electrons plus positrons with a peak around 500 GeV. Finally, there are 
claims that WMAP has observed a 'haze' which could be synchrotron radiation from high energy 
electrons and positrons near the galactic center. This is consistent with dark matter annihilating 
[12] and possibly also with dark matter decaying [13]. 

We will call the decaying dark matter particle ip with mass m^. 

A. Diffuse Gamma-Ray Background From EGRET 

Observations of the diffuse gamma-ray background by EGRET have been used to set hmits on 
particles decaying either into qq or 77 [14], as reproduced in Table 1. We adapted these limits for the 
other decay modes shown in the Tabic. Wc took the age of the universe to be to = 13.72±0.12 Gyr 
4.3 X 10^^ s and the abundance of dark matter to be flBuh^ ~ 0.11 [15]. For the e"'"e~, iJ,'^fJ>~, and 
1/1/ decay modes the limit comes from assuming that these produce a hard W on Z from final state 
radiation ~ 10~^ of the time. The limit on the lifetime is given conservatively as 3 x 10~^ times the 



8 



limit on the W^W~ decay mode, since only one gauge boson is radiated and its energy is slightly 
below For the t+t~ decay mode, the strongest limit comes from considering the hadronic 
branching fraction of the r. The r decays into leptons, eJ/eVr and [xV^Vri 30% of the time. The 
rest of the decay modes have several hadrons and one Vr which carries away at most one-half of the 
energy [16]. Thus we estimate the hadronic fraction of the energy from the decay as | x 0.7 ~ 0.4. 
The limit on qq is relatively insensitive to the mass of the decaying particle in our range of interest 
100 GeV < < 10 TeV. We take this to imply that it depends only on the total energy produced 
in the decay and not as much on the shape of the spectrum, giving a limit on the decay width into 
r's which is a factor of 0.4 of that into qq. To set a limit on decays into the ratio between 

the photon yield from W^W^ and qq is approximated as | from [17]. 



B. Galactic Gamma-Rays From HESS 

HESS observations of gamma rays above 200 GeV from the Galactic ridge [18] can also be used 
to limit the partial decay rates of dark matter with mass > 400 GeV. The limit on the flux of 
gamma rays comes from this HESS analysis in which the flux from an area near the galactic center 
(—0.8° < / < 0.8° and 0.8° <h < 1.5°) was taken as background and subtracted from the flux in the 
galactic center region (—0.8° < / < 0.8° and —0.3° < b < 0.3°) and the resulting flux reported. Our 
limit on the decay mode ^ — 77 is found by taking a similar difference in the flux from decays and 
setting this equal to the observed flux. Because we are considering decays to 77 they give a line in 
the gamma-ray spectrum whose intensity need only be compared to the observed flux in one energy 
bin. This is similar to the analysis in [19] and, as a check, their limit on a dark matter annihilation 
cross section agrees with our quoted limit on the lifetime. 

The photon flux from decays is given by 

riv f 

^decay = -, — - pdrdQ (3) 

where the integral is taken over a line of sight from the earth within a solid angle AQ, r is the 
distance from the earth, is the mass of the dark matter and F is its decay rate, and is the 
number of photons from the decay which we will set equal to 2 for our limits. The density of dark 
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matter is taken as 



Pis) = ^- ^ (4) 

■ 7 / / \ ctN - — - 

U ' 



where s is the radial coordinate from the galactic center. We use the Kravtsov profile [20] with 
(a,/?, 7) = (2,3,0.4), rs = 10 kpc, and po is fixed by p(8.5 kpc) = 0.37 for our limits. This 
gives conservative limits since the flux from the galactic center is much less than in the commonly- 
used NFW profile [21] with (a,/?, 7) = (1,3,1), = 20 kpc, and p(8.5 kpc) = 0.3 There 
is an even more sharply peaked profile, the Moore profile [22], which is defined by Eqn. (4) with 
(a, /3, 7) = (1.5, 3, 1.5), Vg — 28 kpc, and p(8.5 kpc) = 0.27 There is also the very conservative 
Burkert profile [23, 24, 25] 

pM = 7 ^ (8) 



( 



1+f 1 



where po = 0.839 and = 11.7 kpc. We sometimes translate limits on the annihilation cross 



cm- 



section, av, into limits on the decay rate using the fiux from annihilations 



^annihilation 

A conservative limit on decays to 77 is calculated using the Kravtsov profile. If the NFW profile 
is used instead, the limit is stronger r > 5 x 10^^ s \J ^^^y H we had used the Burkert profile the 
hmit would have been much weaker r > 3 x 10^^ s ^J~^^^- In this case the central profile is so 
fiat that there is essentially no difference between the galactic center signal and the fiux from the 
nearby region used for background subtraction, making the limit from our procedure very weak. 
However there would then presumably be a much better limit from just comparing the actual fiux 
at the center (without background subtraction) to what was observed. Thus we believe that the 
quoted limit in Table I is conservative. 

The decay widths to e+e~, jJi^jjr , and t'^t~ can be limited from the HESS observations. These 
light leptons will bremsstrahlung relatively hard photons with a spectrum that can be estimated as 
(see for example [26]) 

'^"-~"^^ + <^--)^Vog('"'-'\-">)r„ (7) 



dx n \ X J \ '^i 
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where Tu is the decay width of ijj — > II, Tuj is the decay width of — > a is the fine structure 
constant, x — — -, and m; is the mass of lepton I. Ignoring the logarthmic dependence on we 
estimate this as giving 10~^ photons per decay for the hght leptons and |10~^ for r's with energy 
high enough to count in the 'edge' feature in the final state radiation spectrum. We then scale the 
limits from decay to 77 by those factors because the edge is assumed to be visible as the line from 
77. Clearly, a more realistic analysis would include a better determination of the observability of 
the edge feature. 

We do not place limits on decays to WW or uu because the spectrum of photons produced 
by final state radiation does not have a large hard component. These can produce many softer 
photons but these arc better limited by lower energy gamma ray observations such as EGRET and 
are counted in the first column. Although the qq decay mode may have a large FSR component, 
this will still give a bound worse than the HESS bound on the 77 mode. Additionally the qq mode 
produces 7r°'s which decay to photons but these are at low energies so HESS cannot place good 
hmits on them. So the qq mode is also better hmited by the EGRET observations. 

C. Neutrino Limits From SuperK, AMANDA, and Frejus 

To find the limits on dark matter decays from astrophysical neutrino observations we start by 
finding the limits on decays directly into two neutrinos using [27]. The given limit on annihilation 
cross section into 1/1/ is almost independent of the dark matter mass in the range 10 GeV < 

< 1 TeV so we simply take it to be a constant. This can be converted to a hmit on the 
decay rate by comparing the neutrino flux from annihilations at the limit to the neutrino flux from 
decays using Eqns. (3) and (6). This gives the limit on the decay rate to vv. This limit comes 
from considering the signal from a 30° half- angle cone around the galactic center (so called "Halo 
Angular" in [27]). The limit from using the signal from the full sky is only a factor of 3 worse. 
Though the signal comes from the direction of the galactic center, the cone is wide enough that the 
given hmit is essentially independent of whether the Kravtsov, NEW, or Burkert proflle is used. 

The limit on the decay rate to e'^e~ is set by the minimum branching ratio for the electrons to 
bremsstrahlung a hard W or Z, which is ~ 10~^. So we conservatively take the hmit to be 3 x 10~^ 
times the limit from neutrino observations on decays into W~^W~ because only one gauge boson 
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is radiated and its energy is slightly less than The limits on jJi^ jJT come from the fact that 
muons always decay as ^ eve^n ^-nd the i/^ tends to carry away almost half the energy. Further 
all produced neutrinos will oscillate a large number of times over these galactic distances before 
reaching the Earth so all neutrinos count equally for detection (and [27] already assumed for their 
limit that all three neutrino species are produced equally). So we ignore the soft V^. and give the 
limit as \ti,i, where Ti,i, is the limit on the lifetime into vv. The factor | comes from assuming 
the has half the energy of the ji and the neutrino background scales as E'"^. Equivalently, the 
bound on Tj,,^ scales with (because the bound on annihilations from [27] is constant in m^) so a 
decay to is like a decay to vv but with a decaying particle of half the mass and half the dark 

matter density. Similarly, the limit on r+r^ of is set by assuming that the Ur, which is always 
produced in r decay, generally carries away about half the energy of the r. Most r decays are two- 
or three-body so we expect this to be a good approximation [16]. 

The W decays one-third of the time to Iv [16] and the neutrino has about half the energy so 
the decay to W^W' is limited to j2'^vv Similarly the Z decays 20% of the time to vv [16] and 
since there are two neutrinos, each of which carries away about half the energy, the limit on decays 
to ZZ is slightly stronger: ^\t^u- 

We do not limit decays to 77 or qq because we expect these to be better limited by direct 
gamma ray and antiproton observations. 

D. Positrons and Antiprotons from PAMELA 

We translate recently pubhshed limits from PAMELA on the annihilation cross sections of dark 
matter into the various final states into limits on the decay rate. This can be done because a dark 
matter particle of mass decaying into one of the given final states (e.g. qq) yields exactly the 
same spectrum of products as two dark matter particles of mass |m^, annihilating into the same 
final state. To translate the limit on annihilation cross section we set $decay = ^annihilation from 
Eqns (3) and (6) but we must use instead of my, in Eqn. (6) for ^annihilation- Also, we integrate 
over the entire sky, AQ = Att, but only over a local sphere out to a radius rmax = 5 kpc. This is a 
crude model for the fact that antiprotons and positrons do not propagate simply like gamma rays 
do. We can ignore the subtleties of this propagation because we are not computing the actual fiux 
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observed, just the ratio between the flux from decays and from annihilations. We then just take the 
simple model that these particles only arrive at earth from a distance of Tmax ~ C(5 kpc). Using 
these assumptions, the lifetime of a decaying dark matter particle with mass that corresponds 
to the annihilation rate of a dark matter particle with mass ^m,p and cross section av is: 

r = 4xl0-s(^)^— ^J. (8) 

It turns out that this is almost independent of halo profile and the size, rmax, of the local sphere 
used to define it. We expect limits from ATIC to be similar to our limits from PAMELA because 
roughly the same signal that fits ATIC will fit PAMELA (see for example Section VI). 

The limits from positrons in the e+e~, n'^ n~ i t~^t~ , and WW channels are translated from the 
annihilation cross section limits from [28]. Really we use the largest annihilation cross section which 
could explain the observed PAMELA positron excess [9] given propagation uncertainties (model B 
of [28]) and translate this into a decay rate. This means that hfetimes around and up to an order of 
magnitude greater than those given in Table I are the best fit lifetimes for explaining the PAMELA 
positron excess. This is most true for the lepton channels, while decays to WW do not seem to fit 
the shape of the positron spectrum very well (see e.g. [17, 28]). Note that these lifetimes are in 
qualitative agreement with those found in [29] . Note that we do not place a limit on the qq channel 
because this is better limited by the PAMELA antiproton measurement. The limit on i/u comes 
from assuming the usual 3 x 10~^ factor times the WW limit from one of the neutrinos producing a 
W from final state radiation. It is possible that a stronger hmit could be set by considering soft W 
bremsstrahlung from the neutrino, turning the neutrino into a hard positron and the limit would 
then come from that positron and not from the decay of the W. We do not attempt to estimate 
this. 

The limit from positrons on the 77 decay channel arises when positrons are produced through 
an off-shell photon. The relative branching ratio of this decay is well known from 7r° decay [30] 

r^^^y Sn \ \mf ) 4/ 

where / is a fermion of charge Q, lighter than ip. Counting the production of muons as well (since 
they always decay to electrons), a conservative estimate for this branching fraction is 4% in our 
range of masses 100 GeV < < 1 TeV. We assume the produced positron to have an energy 
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around ~ The scaling of the hmit on the e^e~ channel with comes from the scaling of the 
number density of dark matter. Thus the energy of the produced positron is not very relevant in 
our range of so the limit on the 77 decay channel is just 4% of the limit on e^e" . 

The limits from antiprotons on the qq and WW channels come from comparing the antiproton 
fluxes computed in [17] with data from PAMELA [31, 32], finding the limit this gives on annihilation 
cross section, and converting this into a limit on the decay rate using Eqn. (8). These are almost 
exactly the same limits as would be derived by translating the limits on annihilation cross section 
from [33] into limits on decay rate. The limits on the e^e~ , yu''"//", t'^t~ , and uv channels come 
from hard W bremsstrahlung from the leptons producing antiprotons and we take this to be 3 x 10~^ 
of the limit on WW (the same factor as used above). Note that these antiproton limits arc most 
applicable in the range 100 GeV < < 1 TeV and are essentially independent of mass in that 
range (equivalently the annihilation cross section limits scale linearly with m^). 

The limit from antiprotons on the 77 decay channel comes when one of the photons is off-shell 
and produces a qq pair. Using Eqn. (9) and summing the contributions from the four light quarks 
(using their current masses) with a factor of 3 for number of colors, we estimate a branching ratio 
B{ip ^qq) ~ 4% at = 100 GeV (rising to ~ 6% at = 1 TeV). Note that this is in rough 
agreement with the one-photon rate found in [34]. If wc conservatively assume that the quark pair 
has energy = then we find that the limit is 0.02 of the limit on the qq decay channel. 

1. Explaining the Electron/Positron Excess 

In the previous section we set limits on the decay rate of a dark matter particle using several 
experiments including PAMELA. The limits from positron observations were less stringent than they 
would have been had PAMELA not seen an excess. In this section we consider what is necessary 
to explain the PAMELA/ATIC positron excesses. A more detailed analysis of individual models is 
presented in Section VL 

Due to the large positron signal and hard spectrum detected by PAMELA the best fit to the 
data is achieved with a decaying (or annihilating) particle with direct channel to leptons. Further, 
the lack of a signal in antiprotons disfavors channels with a large hadronic branching fraction such 
as qq or WW. Thus the PAMELA excess is fit well by a dark matter particle that decays to e+e~. 
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or T+r~ with the hghter two leptons providing the best fit [17, 28]. As in the previous 
section, we translate the cross sections given in [28] using Eqn. (8) to find the range of hfetimes 
which best fit the PAMELA positron data given the propagation uncertainties. We find the best 
fit for masses in the range 100 GeV < < 1 TeV is 



The range comes from the uncertainties in the propagation of positrons in our galaxy and we 
estimate it by using the maximum and minimum propagation models from [28] (models B and C). 
It is difficult to fit just the PAMELA positron excess (not even considering antiprotons) with decays 
to WW unless the dark matter is light < 300 GeV or very heavy > 4 TeV [17]. We will not 
consider these mass ranges because we also wish to fit the ATIC spectrum which requires masses 
in the range 1 TeV < < 2 TeV. Note that while this paper was in preparation [35] appeared 
which generally agrees with the results presented here. 

So for explaining the PAMELA/ATIC excesses, we are most interested in a model in which dark 
matter decay releases between 1 TeV and 2 TeV of energy, dominantly decaying to hard leptons 
with a lifetime given in Eqn. (10). Such a model must also avoid the limits on other subdominant 
decay modes to hadrons which are given in Table I. For example a 1.6 TeV dark matter particle 
produces an endpoint of the electron/positron spectrum around 800 GeV, in agreement with ATIC. 
Such a particle must have a lifetime to decay to leptons in the range 10^^ s < r < 5 x 10^^ s and a 
lifetime to decay to qq longer than r > 10^^ s. 

III. DARK MATTER DECAYS BY DIMENSION 6 OPERATORS 

The decays of a particle with dark matter abundance into the standard model are constrained 
by many astrophysical observations. These limits (see section II) on the decay lifetimes are in the 
range 10^^ to 10^^ seconds. Current experiments like PAMELA, ATIC, Fermi etc. probe even longer 
hfetimes. These hfetimes are in the range expected for the decays of a TeV mass particle through 
dimension 6 GUT suppressed operators. In this section, we present a general operator analysis of 
such dimension 6 operators. 

A dimension 6 operator generated by integrating out a particle of mass M scales as M~'^. This 




(10) 
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strong dependence on M implies a wide range of possible lifetimes from small variations in M. 
The decay lifetime is also a strong function of the phase space available for the decay, with the 
lifetime scaling up rapidly with the number of final state particles produced in the decay Since we 
wish to explore a wide class of possible operators, we will consider decays with different numbers of 
particles in their final state. Motivated by the decay lifetime ~ 10^^ s being probed by experiments, 
we exploit the strong dependence of the lifetime on the scale M to appropriately lower M to counter 
the suppression from multi body phase space factors and yield a lifetime ~ 10^^ s. In Table II, we 
present the scale M required to yield this lifetime for scenarios with different numbers of final state 
particles. The scale M varies from the putative scale where the gauge couplings meet ~ 2 x 10^^ 
GeV for a two body decay to the right handed neutrino mass scale ~ lO^^ GeV for a five body decay. 
We note that the scale lO^^ GeV also emerges as the KK scale in the Horava-Witten scenario. In 
the rest of the paper, we will loosely refer to these scales as Mqut- 



Number of Final 


Scale M (GeV) 


State Particles 




2 


10^6 


3 


3 X 10^^ 


4 


5 X 1014 


5 


1014 



TABLE II: A rough estimate of the scale M that suppresses the dimension 6 operator mediating the decay 
of a TeV mass particle in order to get a lifetime ~ 10^^ s for decays with various numbers of particles in 
the final state. Phase space is accounted for approximately using [36]. Lifetimes scale as M^. Specific 
decays may have other suppression or enhancement factors as discussed in the text. 

The observations of PAMELA/ATIC can be explained through the decays of a TeV mass particle 
with dark matter abundance if its hfetime ~ 10^^ s (see section II). PAMELA, in particular, observes 
an excess in the lepton channel and constrains the hadronic channels. In our operator analysis, we 
highlight operators that can fit the PAMELA/ATIC data. However, wc also include operators that 
dominantly produce other final states like photons, neutrinos and hadrons in our survey. While these 
operators will not explain the PAMELA data, they provide new signals for upcoming experiments 
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like Fermi. For concreteness, we consider SU (5) GUT models. We classify the dimension 6 operators 
into two categories: R-parity conserving operators and R-breaking ones. In the R-conserving case, 
we will add singlet superfield(s) to the MSSM and consider decays from the MSSM to the singlet 
sector and vice versa. The singlets may be representatives of a more complicated sector (see section 
V). In the R-breaking part we will consider the decay of the MSSM LSP into standard model 
particles. As a preview of our results, a partial list of operators and their associated final states are 
summarized in Table III. 



Operator in SU(5) 


Operator in MSSM 


Final State 


Lifetime (sec) 
{Mgut ~ 10^6 GeV) 


Mass Scale (GeV) 
(lifetime ~ 1026 sec) 


R-parity conserving 










stsiotio 


StSQtQ, 5t5C/tC/, S'^SE^E 


leptons 


1026 


10^6 






quarks 


1026 


1016 


Sn0f5\l0f 


S^QL^U, S^UD^E, S^QD^Q 


quarks and leptons 


5 X 1028 


1015 


S%HllOf 


S^LHlE, S^DHlQ 


leptons 


1026 


1016 




S^WemWem, S^WzWz 


7 (line) 


1026 


1016 


Hard R violating 










5/(S5/)5/(I]5/)5/ 


DDDLL 


quarks and leptons 


1037 


1013 


Soft R violating 












4 ~ 
™SUSY ZJ 1) 

'^'^GUT 


quarks 
leptons 
7 + 1/ 


4 X 10^0 
6 X 10^2 
2 X 10^2 


7 X 10" 
1014 
10" 



TABLE III: A partial list of dimension 6 GUT suppressed decay operators. For each operator, we list 
its most probable MSSM final state. The lifetime column gives the shortest lifetime that this operator 
can yield when the scale suppressing the operator is ~ IQi^ GeV. In the mass scale column, we list the 
highest possible scale that can suppress the operator in order for it to yield a lifetime ~ 10^^ seconds. 
Assumptions (see text) about the low energy MSSM spectrum were made in order to derive these results. 
All the operators are in superfield notation except for the soft R violating operators. 
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A. i2-p£irity Conserving Operators 

The lifetimes of the decays of MSSM particles to the MSSM LSP are far shorter than the 
dimension 6 GUT suppressed lifetime ~ 10^^ s currently probed by experiments. In a R-parity 
conserving theory, the MSSM cannot lead to decays with such long lifetimes since the MSSM LSP, 
being the lightest particle that carries R parity, is stable [2]. These decays require the introduction 
of new, TeV scale multiplets. Limits from dark matter direct detection [37, 38] and heavy element 
searches [39] greatly constrain the possible standard model representations of the new particle 
species. In this paper, we will add a new singlet chiral superfield S in addition to the MSSM [77]. 
The singlet may emerge naturally as one of the light moduli of string theory or it could be a 
representative of another sector of the theory (see section V) . 

Decays between the singlet sector and the MSSM can happen through the dimension 6 GUT 
suppressed operators in Table 111 only if there are no other faster decay modes between the two 
sectors. Such decay modes might be allowed if there are lower dimensional operators between the 
singlet sector and the MSSM. In subsection IIIAl, we consider models where lower dimensional 
operators are forbidden by the imposition of a singlet parity under which S has parity -1. This 
parity could be softly broken if the scalar component of the singlet s develops a TeV scale vev 
(5). In subsection III A 2, we consider models without singlet parity that require additional model 
building to ensure the absence of dangerous lower dimensional operators between the singlets and 
the MSSM. 

For the rest of the paper, we will adopt the following SU (5) conventions: S will refer to a SU (5) 
singlet, (5, 5) to a fundamental and an antifundamental of SU{5) and (10, 10) to the antisymmetric 
tensor of SU{5). The subscript / identifies standard model fields, Wa denotes standard model 
gauge fields and Hu and Ha are standard model higgs fields. The subscript GUT refers to a field 
with a GUT scale mass. /, e, z>, q, it and d refer to sleptons and squarks. and refer to 
higgsinos and W refers to a wino. 

1. Models with Singlet Parity 

The R parity conserving operators in Table III that also conserve singlet parity are 
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(11) 



The decay topologies of these operators are determined by the low energy spectrum of the theory. 
SUSY breaking will split and rris, the singlet fermion and scalar masses respectively. It is 
conceivable that SUSY breaking soft masses may make m~ negative leading to a TeV scale vev (s) 
for s. In this case, additional interactions in the singlet sector will be required to stabilize the vev at 
the TeV scale. A decay mode with a singlet vev will typically dominate over modes without a vev 
since the former have fewer particles in their final state leading to smaller phase space suppression 
(see discussion in sub section III Ala). 

Upon SUSY breaking, there are two distinct decay topologies: 

1. A component of the singlet is heavier than the MSSM LSP and this component can then decay 
to it. The relic abundance of the singlet can be generated if the singlet is a part of a more 
complicated sector, for example, through the decays of heavier standard model multiplets (see 
section V) . 

2. Alternatively, the MSSM LSP is heavier than the singlets. In this case, the MSSM LSP will 
decay to the singlets. 

We now divide our discussion further based on the final state particles produced in the decay. 
Motivated by PAMELA/ ATIC and Fermi, we will be particularly interested in operators that 
produce leptons and photons. 

a. Leptonic decays 

The R and singlet parity conserving operators in Table III that contain leptonic final 
states are 




and 




(12) 
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lepton^ 

singlet term ion 




^ siepton 



FIG. 1: A singlet fermion decaying to a lepton, siepton pair with a singlet scalar vev (s) insertion. 

These operators can be generated by integrating out a GUT scale U{1)b-l gauge boson under 
which both the MSSM and the S fields are charged. At low energies the first of these operators will 
lead to couplings of the form 

* Ts^^Z and —^d*s^^d (13) 



while the second operator will lead to similar operators involving u, q and e. Here we have 
suppressed all fiavor indices. 

We first consider the case when the singlet scalar develops a TeV scale vev (s). In this case, 
the operators in (13) mediate the decay of the singlets to the MSSM LSP or vice versa. With a vev 
insertion, the operators in (13) yield: 

Ts^^Z and -^d*s''^d (14) 



When the singlet fermion is heavier than the MSSM LSP, these interactions mediate its decay 
into MSSM states. In particular, when the singlet fermion is heavier than a siepton, it can decay 
to a siepton, lepton pair. The component operators of (14) produce the two body final state (see 
figure 1), 

s^i^ (15) 
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I 



FIG. 2: MSSM neutralino decaying to a singlet LSP and an e^e pair. The decay is dominantly into 

leptons because sleptons are typically lighter than squarks. 

with a lifetime of 

26 /I TeV\ Vl TeV\ V Mqut V 
W.~2xl0»(^) (— j sec (16) 

per lepton generation, where Am — rus — mi. 

When the MSSM LSP is heavier than the singlet fermion, it will decay to the singlet fermion 
and a lepton, anti lepton pair through the operators of equation (12). These decays (see Figure 2) 
are of the form 

LSP + + s (17) 
The lifetime for this three body decay is 

where Ri is the ratio of the LSP mass to the slepton mass and we assumed that ^ mg ; if this is 
not the case should be replaced by the available energy in the decay Am = — mg- Similarly, 
when the singlet fermion is heavier than the MSSM LSP but hghter than the slepton, the singlet 
fermion will decay to the MSSM LSP through a 3 body decay mediated by an off-shell slepton. Note 

that the decay lifetime ~ 10^^ s when the mass scale suppressing the decay is Mqut ~ 10^^ GeV. 
In a U{1)b-l UV completion of these operators, this scale is the vev of the broken B — L gauge 
symmetry. The B — L symmetry must be broken slightly below the GUT scale {i.e. the putative 
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FIG. 3: A singlet scalar s decaying to a slepton pair with a singlet scalar vev (s) insertion. 

scale where the gauge couphngs meet) in order for the three body decays mediated by this gauge 
sector to have hfetimes of interest to this paper. 

The decays discussed above can also produce quarks in their final states. The decay rate is 
a strong function of the phase space available for the decay and is hence a strong function of the 
squark and slepton masses. As discussed in section VII B, the hadronic branching fraction of these 
decays can be suppressed if the squarks are slightly heavier than the sleptons. Since squarks are 
generically heavier than sleptons due to RG running, the suppressed hadronic branching fraction 
observed by PAMELA is a generic feature of these operators. An interesting possibility emerges 
when the spectrum allows for the decay of the singlet fermion to an on-shell slepton, lepton pair. 
This decay produces a primary source of monoenergetic hot leptons. However, the subsequent decay 
of the slepton to the MSSM LSP will also produce a lepton whose energy is cut off by the slepton 
and LSP mass difTcrcncc. With two sources of injection, this decay could explain the secondary 
"bump" seen by ATIC in addition to the primary bump (see sections VII B and VI). 

The scalar singlet s can decay when the singlet gets a vev. In the presence of such a vev, s can 
decay to a pair of scalars {i.e. sleptons and squarks) or fermions {i.e. leptons and quarks). The 
decays of s to a pair of sleptons is also mediated by the operators in (12) which in components yield 
terms of the form 




GUT 



and 




(19) 



as well as terms involving q, u, and e. 
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s will now decay directly to two sleptons (see figure 3) or two squarks, with a lifetime 

.,^,.,,~2xl0-(— ) (— j [^^j sec (20) 

per generation (and per representation). The hadronic branching fraction of this decay is generically 
suppressed since squarks are expected to be heavier than sleptons (see section VII B). s can also 
decay to a pair of fermions. These decays are mediated by the operators 

f^-r/t^/ and -^rS^d (21) 



which can also be extracted from (12). However, the decays of a scalar to a pair of fermions (of 
mass m/) are helicity suppressed by j [78]. These decays will predominantly produce the most 
massive fermion pair that is allowed by phase space. 

Owing to this helicity suppression, the hadronic branching fraction from the decays of s can 
be smaller than 0.1 and accommodate the PAMELA anti-proton constraint if one of the following 
conditions are satisfied by the scalar singlet mass ms, the slepton mass and the top quark mass 
mt. 

• rris ^ nil, ^ ^t- In this case, the spectrum allows both sleptons and the top quark to be 
produced on shell. Hadrons are produced in this process from the decays of the top quark. 
The branching fraction for top production is {j^^ which is smaller than 0.1 if mg ^ Zrrit. 

• rris > 2m[ and rris < 2mt: In this case, sleptons are produced on shell. Hadrons can be 
produced in this process either through direct production of the b quark or through off-shell 
tops. The branching fractions of these hadronic production channels are smaller than 0.1 
since the direct production of b is suppressed by j and the decays mediated by off-shell 
tops are suppressed by additional phase space factors. 

When the scalar singlet does not get a vev, S parity is conserved. Decays between the singlets 
and the MSSM must either involve the decay of the heavier component of the singlet to its hghter 

partner and the MSSM or the decay of the MSSM to the two singlet components. We first consider 
the case when one of the singlet components is heavier than the MSSM LSP. Without loss of 
generality, we assume this component to be the scalar singlet s. s can decay to its fermionic partner 
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and a lepton-slepton pair through a 3 body decay mediated by the operators in (13) if the slepton 
is hght enough to permit the decay. The hfetime for this decay mode is 



/lOOO GeV^ 



4 



The 3 body decay of s to the singlet fermion and slepton-lepton pair may be kinematically 
forbidden if the slepton is heavy, s then decays to the singlet fermion and the MSSM through a 
four body decay: 

s^l+ + l- + S + LSP . (23) 

The lifetime in this case is 

.^flTeVVf Mgut V f Ri^ ^ 



If the MSSM ncutralino is heavier than the singlets, then its decays through the operators in 
(13) are also four body decays similar to the decay discussed above. The lifetime from this decay 
is ~ 10^^ s when Mqut ~ 3 x 10^^ GeV, which is roughly the scale of the right-handed neutrino in 
a see-saw scenario. In fact, if this decay is mediated by a U{1)b-l gauge boson, the scale Mqut 
that suppresses this decay is the vev that breaks the U{1)b-l gauge symmetry which is roughly 
the mass of the right handed neutrino. 

The decays discussed in this section involve decays between the singlet sector and the MSSM 
LSP. In order for these dimension 6 GUT suppressed operators to be involved in the decays between 
these sectors, it is essential that there are no other faster decay modes available in the model. One 
such mode can be provided by a light gravitino. If the gravitino is the MSSM LSP, then the 
superparticles of the MSSM will rapidly decay to the gravitino with lifetimes ~ . If one 
component of the singlet is heavier than the gravitino, then that component will decay to its 
superpartner and the gravitino with a lifetime ~ ^0- where Am is the phase space available for 
this decay. Since we are interested in the decays of TeV mass particles, this scenario is relevant only 
when the gravitino mass is less than a TeV i.e. F < {W^ GeVf. When F < (10^^ GeV)^ 

the above decays occur with hfetimes ~ 10^ s which are far too rapid. 

Another possibility is for the gravitino to be the MSSM LSP and be heavier than the singlets. In 
this case, the gravitino will decay to the singlet sector with a decay rate ~ {^-^^ yielding a lifetime 
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'10" Gevf ^ ^ 



~ 10 s ^-^^ -p — — j which is also far too rapid. The dimension 6 GUT suppressed operators 

discussed in this section lead to astrophysically interesting decays only when the gravitino is not 
the MSSM LSP and has a mass larger than ~ TeV. This forces the primordial SUSY breaking scale 
F ^ (10''^''^ GeV)^, making the gravitino heavier and forcing it off-shell in the decays mediated by 
it. Integrating out the gravitino, operators of the form (^) (^*'^^/^/) generated. The 

masses mg in this operator are the singlet and soft SUSY breaking masses ~ 1 TeV. The decays 
mediated by this operator have lifetimes ^ 10^^ s for F ^ (10^^ GeV)^ and will not compete with 
dimension 6 GUT operators discussed in this section. 

The constraints on the gravitino mass can be evaded if there are more singlets in the theory. 
For example, if the dark sector contains flavor, standard model particles may be emitted during 
"flavor-changing" decays in the dark sector. Consider for example an SU (6) extension of the GUT 
group on an orbifold. The chiral matter flelds and their decomposition to SU{5) representations 
are 

6i = 5i + Si 15i = lOi + 5i . (25) 

where the index i represents flavor. We break the SU{6) symmetry by projecting out the light 
modes of the 5-plet 5 by orbifold boundary conditions at the GUT scale. In the absence of the 
5-plet, the singlet superfields S do not have yukawa couplings that connect them to the MSSM 
fields, eliminating direct decay modes between the singlets and the MSSM. 

Once the heavy off-diagonal SU (6) gauge multiplets are integrated out, operators of the form 

T^SlS,L]L, or j^SlS,D]D, (26) 

are generated. These operators may lead to dark-flavor changing decays of Si Sj + h + Ij or 
Si Sj + di + dj if there are mass splittings amongst the singlets. These splittings can arise due to 
explicit SU (6) breaking terms (which may be present on the brane which breaks this symmetry) . 
Soft SUSY breaking will also contribute to mass splittings in the singlet scalar sector. These 
splittings can cause the decay of a scalar singlet Si to a singlet fermion sj and a lepton Zj, slepton 
Ij pair. In this case, the lepton and slepton emitted in this process will belong to different families. 
The hadronic branching fraction of these operators relative to the leptonic channel depends strongly 
on the masses of the various broken SU (6) gauge bosons. This branching fraction is suppressed if 
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the SU{6) gauge boson that connects the singlets and the leptons is hghter than the boson that 
connects the singlets and the quarks. 

The decays mediated by these operators are immune to the effects of the gravitino since these 
decays explicitly require off-diagonal gauge bosons. A relic abundance of the singlets can again be 
generated through non-thermal processes as discussed in section V. 

b. Decays to Higgses 
The R and singlet parity conserving operators in Table III that contain higgs final states are 



(27) 



These operators are very similar to the operators discussed in subsection III A 1 a. They can also be 
generated by integrating out a GUT scale U{1)b-l gauge sector and the topologies of the decays 
mediated by these operators are also similar to the decay topologies of the operators discussed in 
that subsection. However, since these operators involve final state higgses, they will always yield 
an 0(1) hadronic branching fraction. 

In any particular UV completion, the Icptonic operators discussed in subsection III A 1 a and the 
higgs operators presented in this section may be simultaneously present. The hadronic branching 
fraction of the decays between the singlets and the MSSM is a strong function of the phase space 
available for the various decay modes. If the sleptons are lighter than the squarks and the higgsinos, 
the decays will predominantly proceed via the leptonic channels and hence these UV completions 
will also be compatible with the constraints on the hadronic channel imposed by PAMELA. 

c. Decays to Gauge Bosons 
The only operator in Table III that contains gauge boson final states is 

~IP — ^ ' 

This operator may be generated by integrating out a heavy axion-like or dilaton field that couples 
linearly to both and to S'^. For example, at the GUT scale we may write a superpotential 

W — MgutIOgutIOgut + -^gut-^gut-'^gut + -'^gutIOgutIOgut + -'^gut-S'^ ■ (29) 

where Mqut is a GUT scale mass. Integrating out the heavy IOquts will lead to a one loop coupling 
of the form XouTW^Wa/MGUT (see subsection IV C). Using this effective operator and integrating 
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out Xq\jt leads to the operator 

The decay topology of this operator is similar to the topologies already discussed in subsection 
III A 1 a but leads to new final states. For example, when the scalar singlet s develops a TeV scale 
vev, this operator can lead to decays between the MSSM LSP and the singlet fermion that result 
in the direct production of a monochromatic photon. The lifetime for this decay is 

o .n29^ (5) \V Am \V3 X lO^^GeVX^ 

with Am — nis — ttilsp- 

This signal is noteworthy since Fermi will have an enhanced sensitivity to a photon line up to a 
TeV. Direct decays to monochromatic photons is also a qualitatively different feature permitted for 
decaying dark matter. The production of monochromatic photons from dark matter annihilations 
is loop suppressed and hence annihilations always lead to bigger signals in other standard model 
channels before yielding signals in the photon channel. However, direct decays of dark matter to 
monochromatic photons can happen independently of its decays to other channels. 

The operator discussed in this section will also induce decays involving Zs or decays to a 
chargino and a W boson. The relative rates compared to the photon decay will be set by the 
decay topology, the bino vs. wino composition of dark matter, as well as the relative size of the 
U{1)y vs. SU{2)l couplings. Hadronic decays to a gluino and a gluon are also possible if they are 
kinematically allowed. 



2. S-Numher Violating Operators 

A different class of operators are Kahler terms that break the global S number, for example 

SH^fHilf SHQfH\lQf -^tlO^S^lO; (32) 

The first two operators involve an R-even S and allow the lightest neutralino of the MSSM to decay 
to the fermion component of S. If the neutralino has a significant Higgsino component the decay 
proceeds just through the single dimension 6 vertex and so is three-body. In this case the first 
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operator produces both leptons and hadrons and the second produces only hadrons. In this case, 
it is possible to make the first operator produce mostly leptons if SU(5) breaking effects in the UV 
completion make the coefficient of the S^Hl^LE component dominate over that of the S^HlQD. If 
the neutralino is mostly gaugino the decay proceeds through an off-shell Higgsino, squark, or slepton 
and is four-body. In this case, the first operator could produce both hadrons and leptons. Here, the 
lepton-only channel will dominate if a slepton is lighter than all the Higgsinos and squarks. Again, 
the second produces only hadrons, as in Table III. 

The third operator in Eqn. (32) has an R-odd S, so the neutralino will decay to the scalar 
component of S. This generally produces both leptons and hadrons 

+ + 2 jets Nq^ s + v + 2 jets or A^o ^ « + 3 jets (33) 

If the right handed sleptons are the hghtest sleptons the channel with charged leptons can dominate. 

The decays from the MSSM to S mediated by the operators in Eqn. (32) are either three or 
four body and their rates are as shown in Table III, probing scales as shown in Table II. 

Of course, these operators also allow the S to decay to MSSM particles, if S has a primordial 
abundance. A primordial abundance of S could be generated, for example, by decays at around 
1000 s by dimension 5 operators, as happens in the model in Section VB. The second and third 
operators in Eqn. (32) will always produce hadrons in the decay of but the first operator could 
produce mostly leptons if, for example, the lepton component dominates due to SU(5) breaking 
effects as described above. In a theory with this operator 10/if^5/, if the LSP of the MSSM has 
a large Higgsino component then the decay of the fermion S will be three-body to a Higgsino and 
two leptons. In such a scenario, the scalar s could have two-body decays to lepton or slepton pairs. 
Also possible are three-body decays to a Higgs and two leptons or to a Higgsino, a slepton, and a 
lepton with the slepton then decaying to the LSP and a lepton. Even if the s is lighter than the 
LSP it will still decay just to a Higgs and two leptons or just to two leptons if the channel where 
the Higgs goes to its VEV dominates. The details of the decay depend on the mass spectrum of the 
MSSM particles. Such a decay could yield an interesting electron/positron spectrum. Cosmic ray 
observations could then give important evidence for the mass spectrum of the MSSM, as discussed 
in Section VI. 

One example of a UV theory which can generate the operator lO/i^^lO/ is shown in Fig. 4. 
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FIG. 4: A way to generate the operator S*^ 

We have added two new fields at the GUT scale, IOgut and 24gut and their conjugate fields, IOgut 
and 24qut) order to give them vector- like GUT masses. The superpotential is taken to be 

W = 524gut24gut + 10 /24gutT0gut + i/dlOGUTlOGUT (34) 

This preserves a 'heavy parity' under which the GUT scale fields IOgut and 24gut and their 
conjugates are odd and everything else is even. This ensures no mixing happens between the 
new GUT scale fields and the MSSM fields. Further this preserves a PQ symmetry with charges 
QiH^) = Q{Hd) = Q{24gvt) = 2, Q{lOf) = Q{5f) = g(TOGUT) = -1 and Q{S) = -4. These 
symmetries and R-parity forbid all dangerous operators of dimension lower than 6 that would cause 
a faster decay, except for SHuHd which is a superpotential term and so will not be generated if it 
does not exist at tree level. 

Similar box diagram ways exist to generate the other S'-number violating operators. A combi- 
nation of PQ symmetry and R parity forbids dangerous operators of dimension 5 or lower for the 
three operators in Eqn. (32) except for SH^Hd in the case of the first two operators and S5fHu 
for the third operator. These are in the superpotential and so will not be generated if not there at 
tree level. 
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Another possible UV model to generate these operators is to expand the GUT gauge group 
beyond SU(5) and integrate out the heavy (GUT scale) gauge bosons. For example, the operator 
yS'nO/SjlO/ may be generated by integrating out an S'O(IO) gauge boson at the GUT scale. If 
this is the case, the field S* is a right handed neutrino and some model building would be required 
to assure the decay of the neutralino does not happen by dimension five operators mediated by 
Yukawa couplings. This may happen if the lightest right handed sneutrino has no Yukawa coupling, 
suggesting one of the neutrinos would be completely massless. Generating one of the other operators 
involving a Higgs would require an even larger group than SO (10). We do not consider such models 
further. 

B. R-parity Breaking Operators 

The minimal extension of the SSM allowing the dark matter to decay without introducing 
any new light particles arises if R-parity is broken. R-parity is a symmetry imposed to forbid 
renormalizable superpotential operators, 

UDD, QDL, LLE, and H^L, (35) 

that would otherwise cause very rapid proton decay. As a byproduct, it stabilizes the Lightest 
Supersymmetric Particle (LSP) which, if neutral, makes an excellent dark matter candidate. Con- 
sequently, dark matter decay may indicate that R-parity is broken. 

In this section, we will connect the smallness of R-breaking to the hierarchy between the weak 
and the GUT scales. If R-parity violating effects are mediated through GUT scale particles, their 
effects can be suppressed by the GUT scale. But the SUSY non-renormalization theorem is not 
enough to protect the theory from dimension 4 or dimension 5 R-breaking operators, which would 
lead to too rapid LSP decay; for example, if R-parity is broken and there is no additional symmetry 
replacing it, kinetic mixings, such as H^^L, are allowed and are not suppressed by the high scale. 

To illustrate this point consider a dimension 6 operator H^H^fbflQf. In the presence of the 
MSSM Yukawa interactions, there is no symmetry that forbids the dimension 5 operator lO/if^i^^, 
and prevents it from being generated in a UV completed theory, as can be seen from the existence 
of diagram 5. Note, that in this example there is no problem at the effective field theory level 
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because the diagram in Fig. 5 cannot generate 10 fHlHd- If one starts with an effective field theory 
involving only dimension 6 operators, the low energy loops will never generate dimension 5 terms 
because Lorentz symmetry leads to cancellation of the linear divergencies. So the diagram in Fig. 
5 by itself gives rise only to dimension 6 operators like V^IO fHlHa, but not to dimension 5 terms. 
Still the existence of such a diagram is a signal that it will be challenging to UV complete such a 
theory without generating the lOfH^Hd term as well. 




Hi 

FIG. 5: VllOfHlHd generated by if„i?d5/5/10/ and the MSSM Yukawas 

In fact, there are also examples of operators for which the problem arises directly at the effective 
field theory level. For example, consider the superpotential dimension 6 operator, jj^ — HuLWaW" . 
By closing the gaugino legs (see Fig. 6) it gives rise to the R-parity breaking HuL term. This is a 
superpotential term, so one may think it is not generated. Indeed, this loop is zero in the SUSY 
limit. However, in the presence of a SUSY breaking gaugino mass mi/2 the loop is non-zero and 
quadratically divergent, so it gives rise to the term mi/2HuL. The presence of such a quadratic 
divergence does not contradict the lore that the quadratic divergencies are cancelled in softly broken 
SUSY, because at the end of the day we obtained a mass term of order the SUSY breaking scale 

These problems lead us to two possible ways of consistently implementing R-parity breaking 
implying dimension 6 dark matter decays. The first possibility is to replace R-parity by another 
discrete symmetry that forbids the dimension 4 terms in (35) and also dimension 5 operators that 
give rise to the dark matter decays, but allows dimension 6 decays. An alternative proposal is that 
the R-parity is violated by a tiny amount, that is not put in by hand, but related to the SUSY 
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mi/2 




FIG. 6: Gaugino loop generating LHu from jj^ — H^LWaW^ after SUSY breaking 

™GI7T 

breaking scale. Let us illustrate each of these options in more detail. 
1. Hard R-parity Breaking 

R-parity is not the only discrete symmetry that forbids the dangerous lepton and baryon number 
violating operators (35). Alternative discrete symmetries may arise from broken gauge symmetries 
and insure the longevity of the proton. Heavy fields may have couplings that preserve these sym- 
metries but not R-parity [40, 41]. As a result, these symmetries allow for the LSP to decay at the 
non-renormalizable level. One such operator, 

DDDLL 



M2 



(36) 



arises when there is GUT scale antisymmetric representation of SU(5), while the fundamental 
theory obeys a symmetry (see Table IIIB 1). This operator is generated by the combination of 
the couphngs 

W D 10gc/t5 /5 / + IOgc/tIOgc/tS / + SIOgi/tIOgc/t, (37) 

that violates R-parity. S is the adjoint that breaks SU(5) down to the SM gauge group. It is 
essential that it splits the colored and the electroweak parts of I^guti otherwise R-parity violation 
would come through (5/)^ which is zero. 
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Particle 


Z3 charge 


10/ 




5/ 


1 


Hu 


■4tt 

e 3 


Hd 


j27r 

e 3 


l^GUT 


1 



TABLE IV: The charges of the SM fields and the GUT scale fields under a Z3 discrete that substitutes 
R-parity 

The Z3 symmetry also insures that there are no kinetic mixings between Ught and heavy fields 
that introduce rapid DM decay. The LSP decays to 5 SM fermions through a sparticle loop (see Fig. 
7) 

X ^ 3 jets + 2£ (38) 

The rate of decay is of order, 

'''' 

where we took into account the five body phase space and loop suppression. 

This example illustrates two generic features of LSP decays in the presence of symmetries 
replacing R-parity. Typically the allowed operators have a large number of legs, such that the decay 
rates are significantly suppressed by final state phase space. Such operators also tend to involve 
quarks, resulting in order one hadronic branching fractions. 



2. Soft R-parity Breaking 

Even if R-parity is violated only through GUT fields, in the absence of a symmetry, kinetic 
mixings can still generate order one R-parity violating effects. This problem could be avoided if 

R parity is broken only through SUSY breaking effects involving GUT fields. These effects will be 
communicated to the MSSM through the GUT fields resulting in suppressions ~ ( Jmgvt ^ " T^^®' 
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FIG. 7: LSP decay in a theory with a Z3 symmetry substituting R-parity 



for example, two pairs of heavy 5 © 5 and two heavy singlets, 5*1 and The interactions between 
the GUT and MSSM fields are: 

MgutSi + MgutSI + Mgut^guti^guti + Mgut5gut2^gut2- (40) 

The important point is that the R-parities of these heavy fields are not fully defined with this 
superpotential- Si,5guTi have equal R-parity and S2,5gut2 have opposite R-parity. To break 
R-parity one needs two soft terms for heavy fields, for instance: 

'>tT'SusySiS2 and mlusY^GUTi^GUT2 (41) 

Consequently, the R-parity breaking coefficients in the MSSM sector are proportional to the product 
of the two soft masses and are always suppressed by at least M^^^ . Indeed, the loop of heavy 
fields generates the R-breaking i?/i-term 

ivIgUT 

which is effectively dimension 6. 

Even though this scenario works at the spurion level, it is hard to implement in a full theory of 
SUSY breaking. First, we need to sequester the source of R-parity violation from the MSSM fields 
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but not from the GUT fields. In a toy model with one extra dimension, the MSSM and R-breaking 
fields are located on different branes, while the GUT fields are free to propagate in the bulk. Then 
R-parity breaking is communicated to the MSSM fields only through loops of GUT particles. Care 
has to be taken in order to suppress the effects of gravity that propagates everywhere and may 
communicate unsuppressed R-violation to the MSSM. This is ensured if the soft terms for the 
heavy fields are generated by a gauge mediation mechanism. However, as discussed in sub section 
IIIAl, the MSSM LSP can decay through these dimension 6 GUT suppressed operators only if 
the LSP is not the gravitino. This requirement forces the F-term responsible for the MSSM soft 
masses to be much larger than the one responsible for R-breaking. Consequently, we need two very 
different scales of SUSY breaking, one for the MSSM sector and another for the GUT sector. It is 
not clear if such a SUSY breaking mechanism can be successfully embedded into a UV completion. 



IV. THE PRIMORDIAL LITHIUM PROBLEMS AND DIMENSION 5 DECAYS 

Recent observations [42] of the ^Li/H and ^Li/H ratio in metal-poor halo stars suggest a discrep- 
ancy between the standard big bang nucleosynthesis (BBN) and observationally inferred primordial 
light element abundances. As pointed out in [43, 44], the decay of a particle x with a lifetime 
r ~ 100 — 1000 s could explain this anomaly if the energy density Vtyji^ of x ^"^^ the hadronic 
branching fraction Br{H) are such that Q^h^Br{H) ~ 10~^. The required density of x the 
time of its decay is similar to the expected relic density of a particle of mass ~ 100 GeV with 
electroweak interactions. The decay of such a particle can explain the ^Li and ^Li abundances if 
its hfetime r ~ 100 — 1000 s. Previous work [45] has concentrated on obtaining this lifetime within 
the context of the MSSM either through the decays of the gravitino to the LSP or the decays of 
the NLSP to the gravitino depending upon the low energy SUSY spectrum. Due to the small pro- 
duction cross-section of the gravitino, the former scenario is difficult to test at collider experiments 
and either requires non-thermal mechanisms or tuning of the reheat temperature of the Universe to 
furnish the relic abundance of the gravitino. The latter scenario involving the decay of the NLSP 
to the gravitino has been discussed extensively in the literature [49]. The desire to generate a 
dark matter abundance of gravitinos from the decay of the NLSP while simultaneously allowing for 
the required hadronic branching fraction forces the NLSP to be the stau with a mass ~ 1.5 TeV 
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in the generic parameter space of the theory. The LHC reach for such long hved charged staus 
is ~ 500 GeV [46] making it difficult to test this scenario at the LHC. Other previous work has 
suggested that late decays can affect nucleosynthesis [47, 48] . 

In this paper, wc point out that supersymmetric GUT theories provide a natural home to 
another wide class of models that could explain the Lithium anomalies. Dimension 5 operators 
suppressed by the GUT scale are a generic feature of supersymmetric GUT theories. If x decays 
through such an operator, then its decay rate F ~ ( M"^ ^ — ) naturally ~ (100 s) 



In the following, we perform a general operator analysis of the possible dimension 5 GUT 
suppressed operators that can solve the primordial Lithium abundance problem. For concreteness, 
we consider SU{5) grand unification. A sohition to the primordial Lithium problem involving 
dimension 5 GUT suppressed operators requires the introduction of a new TeV scale particle species 
X- We restrict ourselves to models that involve the addition of complete, vector-like SU (5) multiplets 
in order to retain the success of gauge coupling unification in the MSSM. We will only consider 
operators that preserve R-parity, automatically ensuring the existence of a stable dark matter 
candidate. The operators are classified on the basis of the SU (5) representation of x- In each case, 
we discuss the relic abundance and hadronic branching fraction of the decaying species that solves 
the primordial Lithium problem and the corresponding LHC signatures. 

A. Operator Classification 

We illustrate the possible dimension 5 GUT suppressed operators in Table V. The operators are 
classified on the basis of the SU (5) representation of x- The R-parity of x is indicated by subscripts 
e (R-even) and o (R-odd) and chosen to make the operators in Table V R-invariant. We restrict 
ourselves to the cases when x is a singlet, a fundamental (5, 5) or an antisymmetric tensor (10, 10) 
of the SU (5) group. As illustrated in the second and third column of Table V, it is straightforward 
to construct a number of dimension five operators inducing the decay of x or LSP decay, depending 
which one is lighter, in all these cases. The hadronic branching fraction for the corresponding decays 
is typically quite significant, as the operators involve either Higgs fields or quarks. 

Note that fundamental or antisymmetric representations of SU{5) are not good Dark Matter 
candidates; a stable fundamental would imply dirac DM, which is excluded by direct detection 
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searches [37, 38], while an SU{5) antisymmetric does not have a neutral component. Only when 
X is a singlet is the LSP phenomenologically allowed to be heavier than x and decay to it. In this 
case, X naturally inherits the LSP thermal relic abundance. If this singlet is heavier than the LSP, 
one needs to introduce TeV scale interactions that would give x ^ thermal calculable abundance, 
or rely on non-thermal production mechanisms such as tuning of the reheat temperature, in order 
to yield the required x abundance. Due to this difference between x's that do and do not carry 
SM charges, we have only included the quadratic in x operator in the of a singlet x- This operator 
allows for the LSP decay into x's, but doesn't make the lightest component of x unstable. 

A small subtlety in all the cases is that we can also construct relevant and marginal gauge 
invariant operators involving x and the MSSM fields. If present, they would mediate the decay of x 
without GUT scale suppression. These dangerous operators are collected in the last column of Table 
V. In particular, the kinetic mixing terms such as xl^f (here Xo is an R-odd antifundamental GUT 
multiplet; the effect of this term is equivalent to the mixing between Xo and 5/ in the MSSM Yukawa 
interactions) are not protected by SUSY non-renormalization theorems and will be inevitably gen- 
erated unless forbidden by some additional symmetries. However, it is relatively straightforward to 
impose Peccei-Quinn symmetries that either forbid these operators or make them GUT suppressed 
as well. 

Let us illustrate how this works in several concrete examples. If x is an R-even singlet, the only 
dangerous operator is XeHuH^- We can forbid it by imposing a discrete Z2 symmetry, x ~^ ~X- 
With this symmetry, the only allowed dimension five operator is xlH^H^. As a result, x is the dark 
matter, and the lithium problem is solved by the decay of the LSP into pairs of x particles. 

Other supcrpotential dimension 5 terms for an R-even singlet are of the form x^i^ where Wi are 
Yukawa terms present in the MSSM Lagrangian. In this case, we need to use a symmetry other than 
Z2 under which x is neutral. The Higgs fields carry charges Qu and Qd such that Qu + 7^ 0, and 
the charges of the MSSM matter fields are determined by requiring that the MSSM Yukawa terms 
do not violate the symmetry. Then the MSSM /i term is forbidden by this symmetry, however, it 
can be generated if the symmetry is spontaneously broken by the TeV vev of the field S with charge 
{—Qu — Qd) coupled to Higgses through the term SHuHd- The new massless Goldstone boson 
will not appear if the actual symmetry of the action is just a discrete subgroup of this continuous 
symmetry, so that one can make all components of S massive. The lowest dimension operator 
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involving x, S and MSSM Higgses is SxHuHa- When S develops a vev it gives rise to the marginal 
operator xHuHa mediating x decay, whose coefficient is naturally suppressed by the ratio of the 
soft PQ breaking scale, /x, to the scale where the dimension five operator is generated, Mqut- 

It is straightforward to generalize these arguments to other cases as well. For instance, for 
the fundamental x to avoid the kinetic mixings with Higgses or matter fields one can impose a 
discrete symmetry which is a subgroup of the continuous PQ symmetry with the following charge 
assignments for the MSSM fields 

QiOf — 1, QHd — ~Qhu — 2, Q^f — —3 

The charges are chosen in such a way that the MSSM superpotcntial is invariant. If the charge of x is 
equal to — 8 the dimension 5 operator XeHyhfhf is allowed, while all other operators mediating 
X decay are forbidden. In fact, one can check that this is the only possible charge assignment that 
avoids kinetic mixing and allows x to decay without requiring soft breaking of the PQ symmetry. 
There are more possibilities if the Higgs charges are not opposite so that the PQ symmetry is softly 
broken by the /i term. Finally, let us use this example to illustrate that there is no problem to go 
beyond the effective theory analysis and construct renormalizable models generating the required 
operators at the GUT scale. Namely, let us consider the following renormalizable superpotcntial 

W = 10gut5/5/ + IOcvtHuXc + MgutIOgutIOgut • 

Here (IOgut, IOgut) is a pair of R-even GUT scale fields with the PQ charge Qio = 6. By integrating 
out the heavy fields one obtains the dimension five operator XeHybf^f at low energies. 

In addition to decays involving the chiral supermultiplets of the MSSM, singlet xs can also 
have decays involving the gauge supermultiplets Wa through operators of the form x^qW"'. These 
operators can be generated through an axion-like mechanism where x is the goldstone boson of 
some global symmetry broken at the GUT scale. 

B. Relic Abundance 

The energy density VL^h?Br{H) that must be injected into hadrons is plotted against the lifetime 
T of the decaying particle in Figure 8. The ^Li problem can be solved when the lifetime r ~ 
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X SU{5) Rep. 


Superpotential Terms 


Kahler terms 


Soft PQ breaking 


Singlet 


XeWflOfH^, XelOfbfHd, 


XeW\lOf, XeHlHu, 


(mgut) ^eHuHd 




xl,oHuHd, XolO/5/5/, 




/ \ — 

(mgut j XoHybf 










(5,5) 


XeHu5f5f, XeHuHuHd, 


xJlO/10/, 






XolO/10/%, Xo^fHuHd 


XolO}H.a 




(10, 10) 


XelO/10/i?d, Xel0/5/i/„, 


xilO/5^, xi^fbf 


[mL.) ^^5/5/) 




Xo5/5/5/, Xe5f5fHd 


XoHj>\ 





TABLE V: The possible dimension 5 GUT suppressed operators classified on the basis of their generation 
in the superpotential or through soft breaking of PQ symmetry or through kinetic mixing in the Kahler 
potential. The subscript f denotes standard model families, Wa are gauge fields and Hu, Hd are the Higgs 
fields of the MSSM. The R-parity of x is denoted by its subscripts e and o for even and odd parities 
respectively. 

1 

^ i^ vt h'^Br{H) ) * ■ "^^^ operators described in section IV A can cause decays between the MSSM 
and the x sector with these lifetimes. 

1. Electroweak Relics 

Let us first consider how the Lithium problems can be solved by the colorless components 
of fundamental or antisymmetric representations of SU{5). Standard model gauge interactions 
generate a thermal abundance of the electroweak multiplcts in (5,5) and (10, 10). We focus on the 
standard model operators that are extracted from the SU{5) invariant operators in Table V and 
contain the electroweak multiplets (the lepton doublet L in (5, 5) and the right handed positron 
E in (10, 10)) from the x- These operators can be classified into three categories: operators that 
involve quarks or only contain higgses, operators that are purely leptonic and operators that involve 
leptons and higgs doublets. This classification is presented in Table VI. 

Operators in Table V that contain higgs triplets when the xs are electroweak multiplets have not 
been included in Table VI since the higgs triplets are at the GUT scale and cannot cause a dimension 
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FIG. 8: (Color online) This figure from [44] plots the energy density that must be injected into hadrons 
versus the decay lifetime in order to solve the primordial lithium problems. Decays in the red region solve 
the '^Li problem and decays in the green region solve the ^Li problem. 

5 decay of the electroweak x- For example, with x ^ (10, 10), the only gauge invariant operator that 
can be extracted from xBfbfH^ in Table V is E^DDHj where Hj is the color triplet higgs. This 
operator cannot cause a dimension 5 decay of the and is not listed in Table VI. There are also 
operators in Table V from which standard model operators that belong to more than one category in 
Table VI can be extracted. For example, the operator xlO/S/if^ generates an operator containing 
quarks, E^QfOfHa, and an operator containing leptons and higgses, E^EfLfH^. We include this 
operator in both categories in Table VI. A UV completion of this operator involves integrating out 
GUT scale S'f/(5) multiplets. C (1) SU{J)) breaking effects at the GUT scale (like doublet-triplet 
splitting) can result in one of the operators (say, E^Q jD fHd) being suppressed relative to the other 
i^xEfLfHd). ^ 

The relic energy density f^x^^ these electroweak multiplets is ~ 0.1 (^-^^ ■ The decays 
mediated by the operators in Table VI that involve quarks or only contain higgses have O (1) 
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FIG. 9: The decay of a slepton to a lepton and the LSP. 



slepton 




FIG. 10: The decay of a slepton to a lepton, gauge boson and the LSP through an off-shell chargino. 
hadronic branching fractions. These operators can solve the ^Li problem if the x lifetime is ~ 

2 

100 s (^'W'^ ^ (figure 8). However, these decays cannot solve the ^Li problem. A solution to the ^Li 
problem requires a hadronic energy density injection Q^h'^Br{H) ^ 10~^ around a 1000 s. Colhder 
bounds on charged particles imply that > 100 GeV. Consequently, the relic energy density Q^h!^ 
is greater than 10~^. Due to the O (1) hadronic branching fraction, the hadronic energy density 
injected is also greater than 10~^. This injection is too large and over produces ^Li (figure 8). 

The MSSM products of the decays mediated by the purely leptonic operators in Table VI may 
or may not contain sleptons. For example, in the decay of the fermionic component xino of the 
X, the operator x^5/5/ yields a lepton, slepton pair while the operator x^5/ does not produce any 
sleptons. The hadronic branching fraction of decays that involve sleptons depends upon the MSSM 
spectrum. An O (1) hadronic branching fraction will be produced from the decay of the sleptons 
if the SUSY spectrum contains charginos or other neutralinos between the slepton and the LSP. 
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FIG. 11: X decay to a slepton, lepton and higgs. 

These decays can solve the ^Li problem but not the ^Li problem. When this is not the case, the 
slepton predominantly decays to its leptonic partner and the LSP (figure 9) . The slepton decay can 
directly produce hadrons through off-shell charginos (figure 10). But, these decays are suppressed 
by phase space factors and additional gauge couplings. The branching fraction for these processes 
is ~ 10""^. The leptons produced in these decays could be rs. However, even though the r has an 
O (1) branching fraction into hadrons, the hadrons produced in this process are pions. A solution 
to the primordial lithium problem requires the injection of neutrons [43, 44] and hadronic energy 
injected in the form of pions is ineffective in achieving this goal. The dominant hadronic branching 
fraction in these leptonic decays is provided by final state radiation of Z and W bosons off the 
produced lepton doublets. These bosons decay to hadrons with an O (1) branching fraction. Using 
the branching fraction for final state radiation of Z and W bosons from [52], we estimate that the 
relic abundance and hadronic branching fraction from the decays of a 600 GeV - 1 TeV x satisfy 
the constraint Vt-^h?Br{H) ~ 10~^. These decays can solve the ^Li and ^Li problems if the lifetime 
~ 1000 s. 

The hadronic branching fraction of the decays mediated by the operators involving both leptons 
and higgses in Table VI is model dependent. Naively, these operators should have an (9 (1) hadronic 
branching fraction since the higgs decays predominantly to h quarks. However, the higgs operators 
in these fields can be replaced with the higgs vev (/i), resulting in an effective purely leptonic decay 
mode. These leptonic decay modes produce hadrons through final state radiation of Z and W bosons 
with a branching fraction ~ 10""^ as discussed above. The hadronic branching fraction of operators 
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slepton 



FIG. 12: X decay to a slepton and lepton with a higgs vev (h) insertion. 



with both leptons and higgses is the ratio of the decay rate to processes involving the higgs and the 
rate to processes where the higgs is replaced by (h). For example, the operator xHy^f^f can cause 
the X to decay to a higgs, slepton and lepton (figure 11) with a rate Th ~ ( ig^^i^M^ — ) ■ When the 
higgs is replaced by (/i), this operator causes the x to decay to a lepton, slepton pair (figure 12) 
with a rate F; ~ ( sttM'^ — / ■ T'he decay to a higgs directly produces hadrons and the branching 

fraction for this decay mode is ^ = (^) (^)^ ~ lO'^ ( soefcev )^- ^^^^ ~ 500 GeV, this 
hadronic branching fraction is sufficiently small to allow this decay to solve both the ^Li and ^Li 
problems. However, we cannot replace the higgs field by (h) in every such operator in Table VI 
that has leptons and higgses. After electroweak symmetry breaking, the masses of the electrically 
charged Z+ and neutral Z° components of the lepton doublet in (5, 5) are split. The charged fermion 
l^^ is heavier than the neutral fermion by ~ aMz [53, 54] while the masses of the corresponding 
scalar components arc additionally split by ~ cos (2/5) [55]. When the fermion components 
are lighter than the scalars, all the components of rapidly decay to . In this case, the lithium 
problems can be solved with (5^, 5;^)s only through the decays of . With l^' , the operator xl0ji7„ 
does not lead to any SU{3) x U{1)em invariant operators when Hy, is replaced by (h). The only 
decay mode for the Ij^ that is permitted by this operator is a decay to the right handed positron 
and a charged higgs, resulting in an O (1) hadronic branching fraction. As discussed earlier, the 
decays mediated by this operator cannot solve the ^Li problem but can address the ^Li problem. 
The relic abundance of a singlet x is model dependent. A thermal abundance of x can be 
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generated if the x is coupled to new, low energy gauge interactions like a U{1)b-l- It could 
also be produced through the decays of new TeV scale standard model multiplets or through a 
tuning of the reheat temperature of the Universe. The hadronic branching fraction of the singlet x 
operators in Table V is also model dependent. Operators like xlO/10/iJ„ have an O (1) hadronic 
branching fraction since xQ fUfH^, the only standard model operator that can be extracted from it, 
contains quark fields. However, an operator like xlO/S/i^d yields both xQfDfHd and xEfLfHd- A 
UV completion of this operator involves integrating out GUT scale SU{5) multiplets. O (1) SU{5) 
breaking effects at the GUT scale (hke doublet-triplet splitting) can result in the operator xQfDfH^ 
being suppressed relative to xEfLfH^. Due to these tunable model dependences, the decays of 
singlet xs can solve both 'Li and *^Li problems as long as their relic abundance Q^h^ ^ 10^"'. 

Yet another possibility available in the singlet case is that the MSSM LSP is heavier than the 
R-odd component of x- Then the LSP will decay to x the x abundance will be close to the 
dark matter abundance today. In fact, if no other stable particles are added then x itself will be a 
dark matter particle. This scenario shares many similarities with the scenario where the gravitino 
is the lightest R-odd particle, so that the LSP can decay. If there is no additional mechanism 
for generating x (such as the coupling to new low energy gauge interactions like U{1)b-l) the 
decaying MSSM LSP should have rather high relic abundance, Qisph'^ ^ 0.1, depending on the 
mass ratio between the LSP and x- This makes it somewhat challenging to solve ^Li problem. This 
is achievable though if the LSP is a slepton coupled to the singlet through purely leptonic operators. 

It is worth stressing that a generic property of all our models is the presence of several long- 
living particles with somewhat different lifetimes and masses. There are two sources for proliferation 
of different long-living species. The first is related to R-parity. Indeed, as the x-parity x ~^ ~X 
is broken only by dimension five operators, the lightest particle in the x multiplet is always long- 
lived. However, if SUSY breaking mass splitting between the lightest R-odd and R-even particles 
is smaller than the mass of the MSSM LSP, both of these particles are metastable and will decay 
only through x parity violating dimension five operators. Their presence doesn't change much in 
our discussion. Another reason for the existence of several long-lived particles is that we are adding 
new fields in the complete GUT multiplets that contain also colored components. Let's now discuss 
their story. 
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X SU{5) Rep. 


Quark and Higgs operators 


Purely Leptonic 

Operators 


Leptonic operators 

with liiggs doublets 


Singlet 


XlO fbfHd (= xQfDfHd) , 
xlOflOfH^, x^H^Ha, xlO/5/5/, 
XlO^lO/, xWaW^, xhIHu 


Xl0f5fHai= xEfLfHd), 
X%5^5/, xloJlO/ 


(mgut) ^-^■"^/' 
X^\Hd 


(5,5) 


X^10flOf{= PxQfUf), xHuHuHd, 
xWfWfWf 


(aw) ^^5/' 


X^fHuHd, xW\Hu, 
(m^) (^^^- ^^^'^) 


(10, 10) 


xlOfhHu (= E^QfDfHd) , 
xlOflOfHa, X5/5/5/, x^lO/5^ 


X%5f, /X (^) xlO/ 


Xl0/5/i7„ (= E^EfLfHd) , 

(m^) (^^10/' ^^/5/) , 

xHu5\ 



TABLE VL The classifications of dimension 5 operators based on the standard model operators that can be 
extracted from them. These operators are generated using the electroweak multiplets in x- Operators that 
require the higgs triplet fields are ignored. The subscript f denotes standard model families and H^, are 
the Higgs fields of the MSSM. The R-parity of x is chosen in order to make these operators R-invariant. 



2. Colored Relics 

For the fundamental x we have at least one pair of long-lived quarks Xd,dJ where Xd has quantum 
numbers of the right-handed MSSM cZ-quarks (and, as before, if the mass splitting between R-even 
and R-odds x-quarks is small enough, wc have another pair of long-lived colored particles). For 
antisymmetric x we have long-lived vector-like x-quarks both with quantum numbers of the right- 
handed MSSM H-quark Xu,u and left-handed u and d quarks, Xq,q- 

The evolution history in the early Universe is significantly more involved for colored particles 
[58] and the corresponding relic abundances are border-hne to be uncalculable. The point is that 
in general heavy long-lived colored particles experience two epochs of annihilation as the Universe 
expands. First, they suffer from the conventional perturbative annihilations at high temperatures 
before the QCD phase transition. The resulting relic abundance of the colored particles is at the 
level fly,h^ ~ 10~^ (^iy) ■ This is a very interesting number — as follows from Figure 8, if correct it 
would imply that a long-lived colored particle with a mass in the subTeV range solves both Lithium 



45 



problems. 

However, colored particles experience a second stage of annihilations after the QCD phase 
transition, that can significantly reduce the abundance. Indeed, after the QCD phase transition 
they hadronize — get dressed by a soft QCD cloud of the size of order ~ ^c^cd- plausible 
that when two slowly moving hadrons involving heavy colored particles collide, a bound state 
containing two heavy particles forms with geometric cross-section ~ 30 mbarn. This conclusion is 
somewhat counterintuitive — naively, one may think that the soft QCD cloud cannot prevent two 
heavy particles from simply passing by each other without forming a bound object. The argument, 
however, is that the reaction goes into the excited level of the two x system of the size of order 
~ Aq^jj. At low enough temperatures the angular momentum of such a state is close to the typical 
angular momentum of two colliding hadrons, Li ~ (mxTy^^Ag^j-, so that one may satisfy the 
angular momentum conservation law by emission of a few pions. Assuming that the reaction to 
such an excited level is exothermic the geometrical cross-section appears to be a reasonable estimate. 
After the excited state with two x's forms it decays to the ground level and x's annihilate. As a 
result the relic abundance can be reduced to the values below Q^h'^ ~ 10~^ for a TeV mass particles, 
where they don't affect the Lithium abundance. 

Definitely, many of the details of this story arc rather uncertain at the quantitative level and 
this conclusion has to be taken with a grain of salt. We discuss some of the involved uncertainties in 
more detail in section IV D. Interestingly, in the models we are discussing here, one may avoid going 
into the detailed discussion of this complicated process and be rather confident that the residual 
abundance of the colored particles is close to that given by the perturbative calculation. The reason 
is that in order for the above mechanism to operate the two x particles in the bound state should be 
able to annihilate with each other. This is the case for some candidate long-lived colored particles, 
such as gluino in the split SUSY scenario or stop NLSP decaying into gravitino, but not always 
true. 

For instance, for antisymmetric x annihilation is possible in some of the bound states (e.g., 
XuXu)j but not in the others (e.g., XuXq)- Once formed, these bound states go to the Coulombic 
ground state, which is compact and doesn't get converted into other mesons any longer. Through 

weak decays such a meson decays to the energetically preffered neutral ground state, that survives 
until individual x particles decay. In fact, it is likely that an original x-hadron is a baryon, given 
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that the reaction converting x-meson and ordinary proton into x-baryon and pion is exothermic. 
This doesn't change the story much; after two transitions one obtains in this case baryons containing 
three x's, such as XuXqXq- Again, such a baryon will decay to the stable state and will survive till 
individual x particles decay. Similarly, for fundamental x an order one fraction of them ends up 
being in the compact baryonic state XdXdXd which is safe with respect to annihilations. 

We don't attempt here to analyze the above processes at the precision level, but this discus- 
sion implies that the resulting abundance of colored x's while being somewhat reduced from its 
perturbative value is still high enough to solve ^Li problem, or even both Lithium problems. 

C. Supersymmetric Axion 

We already mentioned dimension 5 decays involving the gravitino as one of the solution of the 
Lithium problems, that does not involve new particles beyond those present in the MSSM. There 
is another well motivated particle in the MSSM that may have similar effects - the axino (see, [34] 
for a detailed discussion of axino properties and cosmology). 

The axion is a well-motivated new pseudo-scalar particle. It solves the strong CP-problem, 
and may constitute a fraction, or all, of the dark matter. Axion-like particles are also generic 
in string models. In supersymmetric models, the axion is a part of a chiral supermultiplet S, so 
it comes together with the scalar (saxino) and the fermionic (axino, a) superpartners. Being a 
(pseudo)Goldstone boson, the axion supermultiplet couples to the MSSM fields suppressed by the 
PQ breaking scale /„. The leading interactions are with the gauge fields 



where Cj are model-dependent coefficients of order one. These interactions are often generated at 
the one-loop level, for instance, by integrating out vector-like fields acquiring the mass from the vev 
of S due to interactions like SQQ (KSVZ model [59, 60]). As a result, as compared to the general 
model-independent analysis above these dimension 5 operators contain extra one-loop suppression 
factors. Consequently, to be relevant for the Lithium problem the high-energy scale entering 
here has to be somewhat lower than the GUT scale, /„ ~ 10^"^ GeV. Still, this scale is intriguingly 
close to the GUT scale. 




+ h.c. , 



(42) 
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The relevance of axino for the Lithium problems crucially depends on its mass. We will focus 
our attention on the gravity mediated SUSY breaking, so that there is no light gravitino. Then the 
saxion receives a mass of the same order as other soft masses, ~ Fgusy/^Pi- On the other hand, 
the mass of axino is a hard SUSY breaking term and its value is highly model-dependent. If it is 
generated at the loop level, it is suppressed by at least one extra loop factor, and varies between 
~GeV down to the keV range. It may also be generated at the tree level after SUSY is broken, if 
SUSY breaking triggers some singlet fields to develop vev's giving rise to the axino mass. In this 
way the axino acquires mass of the same order as other soft masses. 

In all other respects, the axino is a particular example of adding an MSSM singlet. If the axino 
is lighter than the LSP, and the MSSM LSP is bino-like, it will decay to axino and photon. The 
corresponding lifetime is 

/ fa \VlTeV ' 



T ~ 



lO^sec 



where Am^^ is the axino-LSP mass difference. The hadronic branching fraction in this case is due 
to decays with virtual photon producing quarks and is at the few percent level, see Eqn. (9). If 
the resulting axino is the only cold dark matter component now these decays may solve the ^Li 
problem, but not the ^Li problem. If the dominant component of the cold dark matter is the axion, 
and the LSP is light, so that its thermal abundance is low, the LSP decay to the axino may solve 
both Lithium problems. 



D. Split SUSY 

In the Split SUSY framework [56], the SUSY breaking scale is not the TeV but an intermediate 
scale up to 10^^ GeV. All the scalars are at that scale except one Higgs that is tuned to be light. 
The gauginos and Higgsinos are protected by chiral symmetries and they are at the TeV scale, 
giving thermal dark matter and allowing for the gauge couplings to still unify at the GUT scale. 
The gluino of split SUSY can only decay through an off-shell squark (see Fig. 13) and its lifetime 
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is set by the SUSY breaking scale: 



1 m^, . 

rgluino 

' 1287r3 m%^,y 

~ (100sec)-(!^)Y^^ii^^V (43) 
When the SUSY breaking scale is 10^ — 10^*^ GeV, and the gluino lifetime is 100 — 1000 sec, it 
becomes an excellent candidate for solving the Lithium problem(see Fig. 8), if its abundance is 
between fi^/i^ ~ 10"^ - 10"^ 




FIG. 13: Gluino decay through an off-shell squark 

Even though the gluino lifetime is well determined by the squark masses, its abundance is 
difficult to know because it involves strong dynamics [57, 58]. At a temperature of ~ ^^j^g 
gluino will thermally freeze-out. Its cross-section is calculable, since QCD is still perturbative at 
those times, and its abundance is roughly given by ^h? ~ 10^'^ ( T'^xev )^" ^^^r the QCD phase 
transition, the gluino gets dressed into colored singlet states, R-hadrons, which have a radius of 
~ Ag^^. The question is if these extended states can bind to form gluinonium states, g — g, with a 
geometric cross-section, tiAq^j^ ~ 30 mbarn, that will lead to a second round of gluino annihilation 
and suppress its abundance by several orders of magnitude. 

Whether this second round of annihilations happens depends on the details of R-hadron spec- 
troscopy. 

For example, even if the lightest state is an R-tt, the reaction: 

R-TT + baryon — > R-baryon -|- vr (44) 
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may well be exothermic due to the lightness of tt and, if it has a geometric cross-section, most of 
the gluina will end up in R-baryon states. Gluino annihilation now depends on the reaction: 

R-baryon + R-baryon — > g-g + 2 baryons. (45) 

If this reaction is exothermic, it may suppress the gluino abundance by many orders of magnitude. If 
this reaction is endothermic for high angular momentum gluinonium states, it may leave a significant 
number of gluinos in R-hadrons and there is no significant reduction in the gluino abundance. These 
uncertainties render the final gluino abundance incalculable. 

However, the case where there is no second round of annihilations after the QCD phase transition 
provides an interesting scenario for the LHC. The gluino mass range that gives the measured ^Li 
abundance is 300 — 7000 GeV [57]. If the gluino is lighter than ~ 2 TeV it will be produced at 
the LHC, stop inside the detector and then decay, allowing its discovery, the measurement of its 
lifetime, and the determination of the SUSY breaking scale [61]. As the gluino is the heaviest of the 
gauginos, this also suggests that all low scale Split particles have a good chance of being discovered 
at the LHC. For a 300 GeV gluino of a 1000 s lifetime, its decay can solve both the ^Li and ^Li 
problems (see Fig. 8) and it will be abundantly produced at the LHC. The potential discovery of a 
gluino with a lifetime 100 — 1000 sec at the LHC will solidify the primordial origin of the discrepancy 
between the measured Li abundances and sBBN. 



V. MODELS FOR LITHIUM AND DECAYING DARK MATTER 

The primordial lithium abundance discrepancies and the observations of PAMELA/ATIC can 
be explained by the decays of a TeV mass particle through dimension 5 and 6 GUT suppressed 
decays (see sections HI and IV). The dimension 6 operators 

in section III allow decays between a singlet sector Sm and the MSSM with lifetimes long enough 
to explain the signals of PAMELA/ATIC. A large class of dimension 5 operators were discussed in 
section IV. In particular, operators of the form 

lOmSmhh Srr^yVaW^ ^^^^ 

Mgut Mqvt 
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allow for the population of the singlet Sm sector through dimension 5 decays from the MSSM or a 
new TeV scale SU{5) vector-like sector, for example 10^. In this section, we will show that both 
of these decays can be naturally embedded into SUSY GUT models. These models can naturally 

solve both the primordial lithium abundance problem and the observations of PAMELA/ATIC. We 
also consider models yielding mono-energetic photons that give qualitatively new signals for Fermi. 



A. 50(10) Model 

In S'O(IO), the MSSM superpotential is: 



Wmssm = A/16/16/10;, + fxlOhlOh (48) 

where 16/ and 10^ are family and higgs multiplets. We add TeV scale multiplets (16^, 16^) and a 
GUT scale IOgut along with the following interactions: 

W = A16^16/10gut + ml6^16^ + MgutIOgutIOgut (49) 

These interaction terms allow for R-parity assignments -1 for IOgut and -|-1 for 16^ and pre- 
serves a m parity under which 16^ and IOgut have odd parity. Integrating out the IOqut field and 
the SO{10) gauge bosons, we generate the dimension 5 and 6 operators: 




9 , , , and / c^^^ "\ ' ' (50) 

Mgut J J \l6ny \ Mgvt / J \ Mj \ ^ > 



16^16^16/16/\ [ .( 1 Wl6^16^10n ^ /• / 16t^l6U6}16/ 



where Mb-l is the vev that breaks the S'O(IO) U{\)b-l gauge symmetry. The R and m parity 
assignments forbid dangerous, lower dimensional Kahler operators like lOQur^lO/i and 16j„16/. The 
dimension 5 operators in (50) connect the components of the 16^ multiplet that are charged under 
the standard model to the singlet component Sm of the IQrn and the MSSM. However, the only 
operator in (50) that allows for two singlet fields Sm to be extracted from 16^ is the dimension 6 
operator ^ ^^"'^^^^/^^/ ^ which yields ^'^^^^^s^^f ^ ^j^g phenomenology of this model is identical 
to that of the SU (5) x C/(l) b-l model discussed below. The decays of the standard model multiplets 
in IQrn to the singlets and the MSSM fields at ~ 1000 s can solve the primordial lithium 
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abundance problems while the decays between the singlets and the MSSM at ~ 10^^ s can reproduce 
the observations of PAMELA/ATIC. 

B. SU{5) X U{1)b-l Model 

We consider a SU{5) x U{1)b-l model, with U{1)b-l broken at the scale Mb-l near the GUT 
scale. The MSSM is represented in this model by the superpotential: 

Wmssm = y^lQflQfH^ + Ajl0/5///d + txH^Ha (51) 

where 10/ and 5/ arc the standard model generations, and Aj arc the yukawa matrices and 
and arc the higgs fields. We add a GUT scale (IOqut, 10gut)j a TeV scale (10^, 10^) and a 
singlet Sjn to this theory with the following additional terms W in the superpotential: 

W = AiIO^^^IOgut + A210gut5/5/ + MgutIOgutIOgut + mlO^lO^ + rUsSmSm (52) 

The mass terms m and are at the TeV scale and Mgut is at the GUT scale. These inter- 
actions allow R parity assignments +1 for 10m, Sm and IOgut- The superpotential also conserves 
a m parity under which Sm and (10m, 10m) have parities -1. Soft SUSY breaking will contribute 
to the scalar masses and lead to mass splittings between the fermion and scalar components. In 
particular, the singlet fermion mass m{ will be different from the singlet scalar mass m^. Integrating 
out the GUT scale field IOqut and the broken U{1)b-l gauge sector, we get the dimension 5 and 
6 operators: 

/^^« (^) • /^'^ (^) - /- (t^) (^) 

Here the Y represent the other chiral multiplets 10m, 10/, 5/, Hy, and Ha in the model. The gauge 
symmetries of the standard model, supersymmetry and R and m parities ensure that the operators 
in (53) are the lowest dimension operators that connect particles carrying m parity {i.e. Sm and 
10m) and the MSSM. 

Consider the phenomenology of this theory when the mass m of the 10m particles arc greater 
than the singlet masses and mg. The 10m are produced with a thermal abundance due to 
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10m 



\ 

FIG. 14: 10^ decaying to Sm and sleptons. 

standard model gauge interactions. They will decay to the singlets Sm and the MSSM particles 
through the dimension 5 operator in (53) with lifetime r ~ 1000 s { ^^^ Y (loi^Gev)^ where Am 
is the mass difference between the 10j„ and the singlet (see figure 14). Following the discussion in 
section IV, the decays of the electroweak and colored multiplets in the 10^ can solve the primordial 
^Li and ^Li abundance problems. The decays of the 10^ generates a relic abundance of the singlets 
s„i and s„j. Since R and m parities are conserved in this model, decays between the singlets and 
the MSSM have to involve the dimension 6 operators in (53) when the 10^ fields are heavier than 
the singlets. These operators are identical to the R-parity conserving operators discussed in section 
III. Following the discussion in that section, the decays mediated by these operators can explain 
the observations of PAMELA/ATIC. 

A decaying particle can explain the observations of ATIC if its mass is ~ 1.2 TeV (see section VI 
and [10]). We first analyze the case where one component of the singlet has a mass ~ 1.2 TeV and 
this component decays to its superpartner and the MSSM. Without loss of generality, we assume 
that the scalar singlet Sm is the decaying particle with mass ms ~ 1.2 TeV. The singlet fermion Sm 
with a mass and the MSSM LSP with mass Mlsp are hght. A rehc abundance of is generated 
in the early universe from the decays of the thermally produced 10^. The relic energy density of 
the 10m is ~ 0.2 (jy^v)^- 10m decays to both Sm and Sm- However, since Sm is lighter than Sm, 
the branching fraction for decays to Sm is higher due to the larger phase space available for the 
decay. When <^ mg, the branching fraction for decays to Sm is ~ (1 - ^) . The relic energy 
density ^Is^h? of ~Sm is - 0.2 (jf^)' (^) (l - . For m ~ 1.5 TeV and m,- - 1.2 TeV, we get 
Qs^/i^ ~ 10~^. The decays of can explain the observations of PAMELA/ATIC if the lifetime 
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FIG. 15: Axino decay to the lighter axino and two photons. 

~ 10^^ s. This lifetime can be obtained from the dimension 6 GUT suppressed operators in (53) if 
the U{1)b-l symmetry is broken shghtly below the GUT scale with a vev Mb-l ~ 3 x 10^^ GeV. 

The other possible decay topology is for the MSSM LSP to decay to the singlets, with the LSP 
mass Mlsp ~ 1-2 TeV. Stringent limits from dark matter direct detection experiments [37] and 
heavy element searches [39] require the thermally generated 10m abundance to decay to the MSSM 
and the singlets. The 10^ can decay to the MSSM and the singlets Sm through the dimension 5 
operator in (53) if the spectrum permits the decay. Since the 10m has R parity +1, its fermionic 
components 10m can decay only if their mass is greater than the LSP mass ~ 1.2 TeV. The scalars 
10m can decay to the standard model and the singlets as long as the scalars are heavier than the 
singlets. However, if the scalars 10m and the singlets are too light, the MSSM LSP will decay to 
10m and the singlet fermion through the dimension 5 operator in (53). This decay occurs with a 
lifetime ~ 1000 s and will not explain the PAMELA observations. This decay mode must be shut 
off in order for the MSSM LSP to decay through the dimension 6 operator in (53). This is achieved 
if the sum of the 10m and singlet masses are larger than the LSP mass. 

C. Supersymmetric Axions 

The interesting example of the setup that provides both dimension 5 and dimension 6 mediated 
decays is a mild generalization of the axino model for Lithium. Namely, let's assume that there 
are two axion-like particles corresponding to different PQ symmetries, so that they are coupled to 



54 



MSSM through different combinations of operators as in (42) (in other words, coefficients Cj are 
different in the two cases) . 

In this case there are two axinos and after the LSP decays the dark matter will be a mixture 

of the two. The heavier of the axinos is unstable, but its decays involve insertion of two dimension 
five operators, so effectively it has dimension 6 suppression. Typically, the fastest decay channel is 
oi ^ 02 + 27 (see Fig. 15), giving 



where mi2 is the axino mass difference. The decays involving other gauge bosons arc suppressed 
by the higher mass of the intermediate gaugino in the diagram in Fig. 15. As discussed in section 
IV C, axino masses are highly model dependent and may be as high as of order msusY- If that is 
the case for one of the axinos, its decays will produce monoenergetic photons potentially observable 
by Fermi. 

Note that at the one loop level the two axions (and axinos) are mixed by the gauge loops (see 
Fig. 16). This mixing is generated at the high scale, so there is no reason for it to be small. However, 
this does not change the conclusion that the decay of one axino to the other is effectively dimension 
6. Indeed, after one diagonalizes axinos kinetic terms and mass matrices, at the dimension 5 level 
the resulting eigenstates have only couplings of the form (42). An intuitive way to understand this 
is to note that axion is a (pseudo)Goldstone boson, so each axion (and axino) carries a factor of 

with it. Consequently processes involving two of them have dimension 6 suppression. 

In principle, one may generalize this story to the broader class of models explaining Lithium 
problems by dimension 5 decays. Namely, one may consider two singlet fields (not necessarily 
axinos), that do not couple directly to each other and decay through the dimension 5 operators 
listed in the Table V. However, generically, there is a danger that the dimension 5 operator involving 
two singlets can be generated and as a result the transition of one singlet to another will be rapid. 
One can avoid this problem for some choices of operators by imposing additional symmetries, and 
for other choices of operators one may check that there are UV divergent diagrams already at the 
level of effective theory that makes this proposal invalid. The advantage of the axino scenario, is 
that the absence of the dangerous dimension 5 operators is built in automatically. 
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FIG. 16: Mixing between two axinos induced by the gauge loop. 

VI. ASTROPHYSICAL SIGNALS 

In this section we consider the astrophysical signals of dark matter decaying through dimension 
6 operators at current and upcoming experiments including PAMELA, ATIC, and Fermi. The main 
astrophysical signals of long-lived particles decaying at around 1000 s from dimension 5 operators are 
changes to the light element abundances as produced during BBN. As well as solving the Lithium 
problems, such a decay could make a prediction for the abundance of ^Be [44]. These decays 
may also generate a significant component of the current dark matter abundance of the Universe 
(see section V), resulting in the production of naturally warm dark matter that may contribute to 
observed erasure of small structure [70]. 

A. Electrons and Positrons and Signals of SUSY 

The PAMELA satellite has reported [9] a rise in the observed positron ratio above the expected 
background starting at about 10 GeV and continuing to 100 GeV where the reported data ends. The 
ATIC balloon [10] experiment has observed an interesting feature in the combined flux of electrons 
and positrons which begins to deviate from a simple power law around 100 GeV, peaks at around 
650 GeV, and descends steeply thereafter. The HESS experiment [62] has added new data above 
700 GeV, clearly observing the steep decline of the spectrum. It is attractive to explain both the 
PAMELA and the ATIC signals as coming from a new population of electrons and positrons which 
are the products of dimension six decay of dark matter suppressed by the GUT scale. As we have 
noted in section II D 1, such a GUT suppressed decay agrees with the ATIC and PAMELA rates, 
both of order 10^® seconds. In this section we will show that the spectral shape observed by both 
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experiments can be explained as well. We will find that decays of dark matter allow for a much 
better fit to ATIC data, as compared with models considered thus far, when dark matter is allowed 
to decay to both standard model particles and superpartners. 

Beyond the rough agreement in rates, additional important information may be gained by 
studying the spectral shape of the ATIC feature as well as the PAMELA rise. These shapes will 
become even clearer as upcoming experiments will add new data. The PAMELA experiment is 
expected to release positron data up to ~ 270 GeV as well as measure the combined electron 
positron fiux up to 2 TeV. The Fermi sateUite, though optimized for gamma rays, has a large 
acceptance and is also capable of measuring the combined fiux. HESS is also capable of measuring 
the combined flux and may be able to add to their already published by measuring the flux at 
energies down to ~ 200 GeV, checking the ATIC bump. As has been shown for annihilations [63], 
assuming a new smooth component of positrons that would explain the PAMELA rise accompanied 
by an equal spectrum of new electrons and extrapolating to higher energies, gives rough qualitative 
agreement with the height and slope of the ATIC feature. Here we would hke to stress that analysis 
of the precise shape of electron-plus-positron fiux, particularly as data improves, may allow for 
surprising discoveries. This may already be demonstrated by examining ATIC data closely and 
considering models that fit its shape. 

In what follows, we will take the published ATIC data with only statistical errors as an example, 
and assume for now that there are no large systematics. Thus, we will discuss in detail how the 
spectrum may be fit, including its various features. Such an assumption may well be wrong. ATIC 
may have many systematic effects in the data, rendering small features meaningless. We will take 
it seriously merely to give an example of all the information that may be extracted from a full 
spectrum, including small features on top of the overall bump. From this we learn an interesting 
new qualitative feature, that dark matter decays can generically produce several features in the 
spectrum and not merely a smooth rise and fall around the dark matter mass. However, we do not 
beheve that the ATIC data are conclusive yet. Future data with better systematics, including from 
the Fermi satellite, are necessary before we can definitively conclude anything from the details of 
the electron spectrum. 

The ATIC spectrum, shown in figure 18, may in fact contain two features - a soft feature around 
100-300 GeV and a hard feature at 300-800 GeV. The hard feature may extend to yet higher energies 
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if HESS data is considered. The collaboration itself and the following literature has generally 
highlighted the hard and more pronounced feature, however it should be noted that the statistical 
errors in the 100-300 GeV range are quite small, and when compared with a power law background 
extrapolated from low energies, the soft feature may well be more significant statistically. The soft 
feature also has a peculiar shape; a sharp rise just below 100 GeV, followed by a nearly flat, or 
slightly descending spectrum up to 300 GeV. The shape of the hard feature is less clear, and may 
be interpreted either as ending sharply at 800 GeV or ending more smoothly around 2 TeV if HESS 
data is added. It is tempting to fit the two ATIC features by two different products of the same 
decay which would give these shapes. 

If dark matter is a singlet which is heavier than the MSSM LSP, the products of its decay may 
be superpartners. If this is the case, the decay of dark matter may potentially allow us to discover 
supersymmetry and measure its spectrum. Consider for example the simple dimension six operator 
of equation (12), which allows the singlet fermion to decay to an electron and a selectron with a 
lifetime given by 

.6 / l TeVV /I TeVV / Mgut 



2x10^6 ^-^ ^TT^ sec. (54) 

This is a two body decay and both particles are mono-energetic. In particular the electron (or 
positron) energy will be 

ml -ml 

The selectron, being a scalar will then decay isotropically in its rest frame. Assuming the decay is 
directly to an electron and a neutralino LSP the isotropic decay will give a fiat energy distribution 
in the "lab" frame between two edges 

msjml-mlsp) ml - mjgp 

^ 2^, ^ 2m. ^^^^ 

The combined spectrum of the electron-positron pair emitted in the decay is shown schematically 
in figure 17. Note that if one were to measure the three energies Eg, E_ and the mass of the 
dark matter, the selectron and the LSP may be solved for without ambiguity. We further notice 
that the fiat spectrum that the selectron produces is reminiscent of the plateau above 100 GeV in 
the ATIC fiux. The hard monochromatic electron may produce the hard ATIC feature. 



58 



6q 
-a 



selection 



electron 



E 



+ 



E, 



Energy 



FIG. 17: A schematic plot of a spectrum of the electron-positron pair emitted when a singlet dark matter 
particle decays to an electron and a selectron. Such a spectrum provides an explanation of the two features 
in the ATIC flux. 

In the top left panel of figure 18 we show that such a decay can indeed fit the ATIC data 
remarkably well. We have used GALPROP (v5.01) [64] to generate electron and positron back- 
ground and also to inject and propagate the signal, modifying the dark matter annihilation package 
to one for decays. In order to reproduce the ATIC shape, including its sharp rises in flux, we 
used a propagation model with a relatively thin diffusion region, "model B" of [28] (L = 1 kpc) 
which was found to agree with cosmic ray data. In addition, to the propagation model the ATIC 
flux depends on the electron background spectrum which is still highly uncertain [65] (a slope of 
~ —3.15 ±0.35), because it can only be constrained by the electron data itself as well as by photons 
to some extent [79]. This is in contrast with the positron background which may be indirectly linked 
to spectra of nuclei. In this case we have taken a background electron flux with a slope of —3.3 at 
20 GeV. 

Considering ATIC data alone (we will discuss HESS later), the spectrum that gave a good fit (by 
eye) was {Ee, E_, E^) = (700, 100,200) GeV which may be inverted to produce a heavy spectrum 
of 

m, ~ 2800 GeV mg ~ 2000 GeV misp ~ 1800 GeV . (57) 
It is remarkable that such a good fit to both ATIC features may be achieved by assuming the same 
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FIG. 18: (Color online) Left: The combined electron-positron flux for a heavy dark matter candidate 
decaying to an electron and a selectron. The selectron subsequently decays to a neutralino and a positron 
giving a two features in the observed flux. The dotted curves are the various components of the flux - 
electron, selectron and background. ATIC (red circles) and HESS (blue squares) data are shown. The 
darker grey band is the HESS systematic error and the lighter grey is the area covered by this error after 
including the uncertainty in the energy scale shift in both directions. Right: The positron fraction for the 
same decay and PAMELA data. The dotted line is the expected background. 



rate of injection for both the hard electron and the selectron, strengthening the case for them to 
come from the same decay. Assuming the singlet s makes up all of dark matter we can deduce its 
lifetime from the rate that fits th ATIC spectrum to be 1.5 x 10^^ seconds. Using equation (54), 
this decay is probing a scale of 0.8 x 10^^, remarkably close to the GUT scale. 

The heavy selectron is required so that the selectron will be mildly boosted producing the 
feature at relatively low energies. The superpartner spectrum is quite sensitive to the values of 
{Ee, Ej^) and thus to the shape of the combined electron-positron flux. The uncertainties on 
these masses are thus large, being affected by propagation uncertainties as well as the statistical 
and systematic errors of ATIC, HESS and future data. It will be interesting to perform a more 



60 



systematic analysis as data improves. 

It is also interesting to consider other ways to get a similar spectrum. For example, if dark 
matter is a scalar singlet, it may be decaying both to an e^e~ pair and to a e'^e~ pair. This would 

be given for example by the operator 5*^ from Section 111 A 2 where the Higgs gets a vcv. 

This operator docs not have a helicity suppression because the scalar s decays to two left-handed 
spinors. In order to reproduce the hard ATIC feature at 700 GeV and the soft feature between 100 
and 200 GeV, the spectrum is lighter 

m~s ~ 1400 GeV mg ~ 660 GeV rriLSP ~ 500 GeV . (58) 

However, if the operators that lead to decays to electrons and to selectrons arc suppressed by the 
same high scale, as would be expected in a supersymmetric theory, one would expect that the decay 
to electrons would be more rapid due to phase space. The ATIC spectrum on the other hand 
requires the two rates to be equal, which would imply the two rates must be tuned independently 
somehow. 

On the top right panel of figure 18 we show the positron fraction this decay, s e^e^, produces. 
Qualitatively the positron fraction agrees with that seen by PAMELA, a monotonic rise. However, 
the positron fraction in this case stays flat until roughly 40 GcV before climbing rather abruptly 
whereas the PAMELA data goes up more smoothly. Though the fit is not perfect, the quantitative 
disagreement is not very significant. At yet higher energies, beyond the range of current data, the 
positron fraction undergoes a wiggle, transitioning from the soft to the hard components of the 
spectrum. This feature is an interesting prediction for future PAMELA data or other experiments 
such as AMS2. 

The mild disagreement between the PAMELA data and the spectrum that fits ATIC exemplifies 
that there is a mild tension between the detailed shape of the spectra observed by the two experi- 
ments. This tension may be phrased model independently by noticing that given a smooth power 
law for the electron background, say the dotted background fine on the right panel of figure 18, the 
difference between it and the ATIC spectrum may be interpreted as signal, half of which is electrons 
and half of which is positrons. This new component of positrons can be combined with the positron 
background and the total flux to give a positron fraction which has the peculiar shape seen in the 
right panel of figure 18. The mild tension is thus between the two data sets, regardless of whether 
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FIG. 19: (Color online) Left: The combined electron-positron flux for a heavy dark matter candidate 
decaying to an electron-positron pair. The dotted curves are the various components of the flux - signal 
and background. ATIC (red circles) and HESS (blue squares) data are shown. The darker grey band is the 
HESS systematic error and the lighter grey is the area covered by this error after including the uncertainty 
in the energy scale shift in both directions. Right: The positron fraction for the same decay and PAMELA 
data. The dotted line is the expected background. In this decay channel the fit to ATIC data is inferior 
to that in the previous case, but the fit to PAMELA data is slightly improved 

the signal originates from annihilation or decays. 

An alternative approach is to begin with a smoother signal that reproduces the PAMELA shape 
very well, but to give up on the various ATIC features and treat ATIC as a single smoother bump. 
This is shown in figure 19 in which both the combined fiux and the positron fraction are shown for 
a dark matter with a mass of 1.4 TeV that is decaying to an electron-positron pair. The lifetime of 
dark matter in this case is 3 x 10^^ seconds which may result by a dimension six operator suppressed 
by a scale of 1.5 x 10^^ GeV. This spectrum does not require sharp features and a more standard 
propagation model was used, (the diffusive convective model of [66]). 

Though the hard ATIC peak by itself fits a sharp electron positron peak at around 700 GeV, 
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when combined with recent HESS data the hard feature may be interpreted as a broader bump, 
extending to a couple of TeV. The ambiguity in interpretation arises from the uncertainty of whether 
HESS is seeing the background beyond the ATIC bump, or the dechne of the bump itself to a 
lower and steeper background. Though the HESS statistical errors are small, these should be 
added to a larger systematic error (shown as a grey band in the figures), and to an overall energy 
scale uncertainty (whose effect is roughly shown in the figure as a wider and lighter grey band). 
Nonetheless, the HESS decline motivates fitting the harder ATIC feature with a smoother spectrum. 
This may be done without losing the good fit to the soft feature given by a mono-energetic selectron. 
In figure 20 we show the combined flux and positron fraction from a heavy dark matter singlet 
decaying to a mono-energetic muon, with an energy of 2 TeV, and a slow slepton, that produces a 
box-like spectrum as in figure 17. The shape of the hard ATIC peak including the HESS decline 
is qualitatively reproduced, within the propagation uncertainties. This particular spectrum may 
be achieved if a very heavy dark matter candidate (of order 8 TeV) is decaying to a heavy slepton 
(at 5.6 TeV) which decays to a shghtly lighter neutrahno (5.5 TeV). Though this spectrum is very 
heavy, is remarkable that a good flt to both features may be gotten from a single decay. 

B. Gamma-Rays 

Gamma rays can provide the best evidence that dark matter is the explanation for the observed 
electron/positron excesses. Gamma rays produced in our galaxy at energies of hundreds of GeV 
are not bent or scattered and so provide clean directional and spectral information. Such high 
energy gamma-rays have been searched for with Atmospheric Cherenkov Telescopes (ACTs) such 
as HESS and will be searched for in the upcoming Fermi telescope (GLAST). Gamma-rays from 
dark matter decays or annihilations will appear as a diffuse background with greater intensity in 
directions closer to the galactic center. There is also an expected astrophysical diffuse background 
and so to distinguish the dark matter signal may require spectral information as well. The classic 
signals of dark matter are lines, edges, or bumps in the gamma-ray spectrum. 

Dark matter annihilations do not generically produce gamma rays as a primary annihilation 
mode. By contrast, dark matter decays may directly produce photons as a primary decay mode, 
as discussed in Section 111. Operators such as S'^WaW^ or "^W'^^ H^WSbi cause decays of 77 
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FIG. 20: (Color online) Left: The combined electron-positron flux for a heavy dark matter candidate 
decaying to a hard muon (2 TeV) and a selectron. The selectron subsequently decays to a neutralino and 
a positron giving a two features in the observed flux. The dotted curves are the various components of the 
flux - muon, selectron and background. The more rounded spectrum of the hard muon is in qualitative 
agreement with HESS data. ATIC (red circles) and HESS (blue squares) data are shown. The darker grey 
band is the HESS systematic error and the lighter grey is the area covered by this error after including the 
uncertainty in the energy scale shift in both directions. Right: The positron fraction for the same decay 
and PAMELA data. The dotted line is the expected background. 

or ip '-fu which produce monoenergetic photons and neutrinos. This will appear as a line in the 
gamma-ray spectrum, easily distinguishable from backgrounds if the rate is fast enough. This would 
also be a line in the diffuse neutrino spectrum, perhaps detectable at upcoming experiments such as 
IceCube, but it would likely be seen first in gamma rays. The exact reach of upcoming experiments 
depends on the energy of the line as well as the halo profile of the dark matter. However, it is clear 
from the conservative limits in Table I that HESS observations are already in the range to detect 
such decays. Fermi and future ACT observations will significantly extend this reach, probing GUT 
scale suppressed dimension 6 operators. 
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FIG. 21: The spectrum of final state radiation from a = 1600 GeV dark matter particle decaying in 
the galaxy with lifetime r = 10^^ s as would be seen at Fermi with a field of view = 1 sr near the galactic 
center. The black curve is for the decay channel — > e"'"e~, the dashed is for decays to fji^ lJi~ ■ Gray is the 
expected background. 

Even if the primary decay mode does not include photons, they will necessarily be produced by 
final state radiation (FSR), i.e. internal bremsstrahlung, from charged particles that are directly 
produced. In particular, if dark matter decay explains the electron/positron excesses then the 
decay must produce a large number of high energy electrons and positrons. These decays will then 
also produce high energy gamma-rays through FSR with a spectrum given in Eqn. (7). Such 
decays could come, for example, from the operators from Section III. Figure 21 shows the spectrum 
of final state radiation expected from decays of dark matter in the galaxy to e'^e" and as 
would be seen by Fermi. We take the effective area = 8000 cm^, the field of view = 1 sr annulus 
around the galactic center but at least than 10° away from it, and the observing time = 3 yr. This 
spectrum exhibits an 'edge' at half the dark matter mass which could allow these gamma-rays to 
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be distinguished spectroscopically from the diffuse astrophysical background. 

The third general mechanism which produces photons from dark matter decays (or annihila- 
tions) is hadronic decay (e.g. to qq) producing tt'^'s which then decay to photons. These generally 
yield a soft spectrum, similar to the expected diffuse astrophysical backgrounds. Thus, we will not 
consider these signals, as a definitive detection could be challenging. However it is worth noting 
that these signals can be quite large because they are not suppressed by the factor of ~ 10~^ that 
suppresses final state radiation. 

Fig. 22 shows an estimate for the sensitivity of the Fermi telescope to dark matter that decays 
either to 77 or e'^e". The sensitivity to decays to e'^e~ comes observing the edge in the gamma-ray 
spectrum from final state radiation, as in Fig. 21. Note though, the plot shows the reach of Fermi 
in the lifetime of the primary decay ip — » e+e^ and not the FSR decay ip — >■ e~^e~'y. Roughly we 
see that the reach for e"'"e~ is a factor of 10^ worse than for 77, which is the expected probability 
to emit final state radiation. The sensitivity to the gamma ray line signal {ip — > 77) was found by 
demanding that the signal (which is entirely within one energy bin) be larger than 3 times the square 
root of the expected diffuse astrophysical background in that bin. The signals were calculated from 
Eqns. (3) and (7). The expected background was estimated as [26, 67, 68] 



The energy resolution of Fermi was estimated from [69] as the given function below 300 GeV and 
30% between 300 GeV and 1 TeV. The sensitivity to the edge feature was estimated (similarly to 
[26]) by requiring that the number of signal photons within 50% of the energy of the edge be larger 
than 5 times the square root of the number of background photons within that same energy band. 
This sensitivity increases as the size of the energy bin used is increased and docs not depend on 
the actual energy resolution of the detector (except for the assumption that it is smaller than the 
energy bin used). This is clearly only a crude approximation of the actual statistical techniques 
which would be used to search for an edge. 

We exclude a 10° half-angle cone around the galactic center from our analysis because the 
diffuse gamma ray background is expected to be much lower off the galactic plane. The annulus we 
consider, though it overlaps the galactic plane, is a good approximation to the dark matter signal 
from off the galactic plane. Ignoring the galactic center reduces the signal. To be conservative 
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FIG. 22: A guess for the sensitivity of Fermi. The soHd Hnes are the expected reach of the Fermi telescope 
in the lifetime to decay into the modes ?/' — > 77 and ■0 e^e~ as a function of the mass of the decaying 
dark matter. The signal from the second mode comes from internal bremsstrahlung gamma-rays from 
the electrons (but the plotted reach is for the primary decay mode into just e+e"). The sensitivities are 
conservatively estimated using the Burkert (isothermal) profile, though the NFW profile gives essentially 
the same result. The grey band is the region of lifetimes to decay to e^e~ that would explain the PAMELA 
positron excess from Eqn. (10). The ATIC excess would be explained in this same band for masses near 
K, 1500 GeV. 

we use Eqn. (59) for the background everywhere, though it is probably an overestimate for the 
background off the galactic plane. Fermi may also be able to resolve point sources of gamma rays, 
further reducing the astrophysical diffuse background. 

From Fig. 22 we see that the most of the range of lifetimes which could explain the PAMELA 
and ATIC excesses as a decay to e+e" should be observable at Fermi from the final state radiation 
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FIG. 23: (Color online) The flux of galactic gamma rays versus the angle from the galactic center for 
dark matter decay to 2 photons with = 1 TeV and r = 4 x 10^^ s versus dark matter annihilation to 2 
photons with = 500 GeV and av = 3 x 10^^^^^. Note that the spike at angle is cut off by a finite 
angular resolution taken to be 3 x 10~^ sr. Solid lines are the fluxes from decays, dashed lines are from 
annihilations. Listed in order from top to bottom on the left edge of the plot: Moore profile (red), NFW 
profile (blue), Kravtsov profile (black), Burkert or isothermal profile (green). 

produced. The same conclusion would hold for decays to ji^ ji~ . However, decays to three body 
final states would soften the produced electron and positron spectrum, thus softening the spectrum 
of FSR gamma-rays. This signal could be more difficult to observe with Fermi, though the injected 
electrons and positrons must always be rather high energy in order to explain ATIC so there is a 
limit to how soft the gamma-ray spectrum could be for any model which explains ATIC. 

It may be possible to distinguish dark matter decays from annihilations, so long as the dark 
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matter-produced component of the total diffuse gamma ray background can be distinguisfied from 
the spectrum. The intensity of dark matter produced gamma rays can then be measured as a 
function of the observing angle. Because the decay rate scales as n while the annihilation rate 
scales as n^, the dependence of the gamma ray intensity on the angle from the galactic center is 
different for decays and annihilations for a given halo profile. Of course, uncertainty in the halo 
profile could make it difficult to distinguish between the two. In Fig. 23 we have plotted the galactic 
gamma ray flux as a function of the angle from the galactic center at which observations are made 
(these halo profiles are spherically symmetric around the galactic center so only the angle from 
the center matters). The intensity is shown for both decays and annihilations to 2 photons from a 
variety of halo profiles. We have chosen to scale the rates to the standard annihilation cross section 
{av = 3 X 10~^^^^) and the decay rate that corresponds to this from Eqn. (8) (r = 4 x 10^^ s). 
The shape of these curves is independent of the overall normalization. The decay curves in Fig. 23 
exhibit a universal behavior for large angles > |, independent of the halo profile. This stems from 
the fact that the integral of n is essentially just the total amount of dark matter and is relatively 
insensitive to the distribution. Note that this universal shape of the decay curve is significantly 
different from the shape of the annihilation curves. With enough statistics, this difference would 
be readily apparent. The annihilation curve from the Burkert profile is most similar to the decay 
curve and so would be the most difficult to distinguish. 

A telescope such as Fermi is ideally suited to the task of measuring the angular distribution of 
the gamma-rays because of its large, 0{1 sr), field of view and coverage of the entire sky. Further, 
astrophysical backgrounds are significantly lessened away from the galactic ridge, giving an advan- 
tage to Fermi over the existing HESS observations which are mostly dominated by backgrounds. 
Fermi has an acceptance of ~ 3 x 10^^ cm^ s sr. Thus we can see from the scale in Fig. 23 that 
at least with a relatively strong decay mode into photons, Fermi could have enough statistics to 
differentiate the decay signal from an annihilation signal. 

If the WMAP haze is due to dark matter decays producing electrons and positrons in the galactic 
center it may be possible to detect inverse Compton scattered photons with Fermi. A study has 
been done for annihilating dark matter [71]. It may be interesting to check both the possibility of 
explaining the original WMAP haze and of observing the corresponding inverse Compton gamma 
rays at Fermi, for decaying dark matter as well. 
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VII. LHC SIGNALS 



In this section, we study the LHC signals of the dimension 5 and 6 GUT suppressed operators 
considered in this paper. The dimension 5 operators lead to particle decays with lifetimes ~ 
100 — 1000 s and these can have striking LHC signatures, particularly when the decaying particle 
is colored or has electric charge. The dimension 6 operators lead to particle decays with lifetimes 
~ 10^^ s and these decays will not be visible at the LHC. However, the requirement that these 
operators explain the observations of PAMELA/ATIC leads to constraints on the low energy SUSY 
spectrum and these constraints can be probed at the LHC. In the following subsections, we analyze 
the signals of these two kinds of operators. 

A. Signals of Dimension 5 Operators 

The dimension 5 operators introduced in section IV mediate decays between the MSSM and a 
new TeV scale x- The x is in a vector-like representation of SU{b). In the following subsections, 
we analyze the cases when x is a (10, 10), a (5,5) or a singlet of SU{5) and summarize our results 
in Table VII. 

1. Decouplet Relic 

The most striking coUider signals arise in the case when x is a (10, 10). In this case, all the 
particles in the multiplet are either colored or carry electric charge. When these particles are 
produced at a collider, they will leave charged tracks as they barrel out of the detector. This 
signal should be easily visible over backgrounds at the LHC [46]. A significant fraction of these 
particles will stop in the calorimeters and their late time (~ 100 — 1000 s) decays will lead to energy 
deposition in the calorimeters with no activity in either the tracker or the muon chamber further 
enhancing the visibility of this signal [61]. The LHC reach for such colored or electrically charged 
stable particles was analyzed in [46] and found to be ~ 1 TeV for a right handed positron and at 
least 2 TeV for colored particles. As discussed in sub section IV B 1, the primordial ^Li problem 
can be solved through the decays of a electroweak multiplet only if its mass Me^ ^ 1 TeV. In this 
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scenario, the colored multiplets in the x could also be hght enough to be produced at the LHC. 

The production of these colored multiplets allows for more dramatic signatures at the LHC. Let 
us first consider the case when the fermionic components of the colored SU (2) doublets {U, D) in 
(10, 10) are lighter than the scalar components, so that the scalar components, if produced, decay 
rapidly to the fermionic components. The masses of the components of the colored SU (2) doublet 
fermions {U^ , D-^) in (10, 10) are split by ~ aMz ~ 350 MeV due to electroweak symmetry breaking 
[53, 54]. When the heavier component of the doublet is produced at the LHC, it will decay to the 
hghter one with a hfetime ~ (350^ Me v)^ ~ 4 x 10~^^ s = 1.2 cm, producing a displaced vertex 
that releases an energy ~ 350 MeV resulting in the production of soft pions, muons and electrons. 
The lighter component will barrel through the detector where it will leave a charged track and may 
eventually stop and decay after 100 — 1000 s. While it is difficult to trigger on the soft particles 
produced from the displaced vertex [53, 54] , the detector will trigger on the charged track produced 
from the decay and this track could potentially be used to identify the displaced vertex and the soft 
particles produced with it. The observation of the displaced vertex, while experimentally difficult, 
will unveil the presence of the SU{2) colored multiplet. 

It is also possible that the scalar components of the colored doublets are lighter than their 
fermionic partners. In this case, the fermionic components will decay rapidly to the scalar com- 
ponents. The contribution from electroweak symmetry breaking to the mass splitting between the 
components of the colored SU (2) scalar doublets is ~ cos (2(3) [55] . The decays of the heavier 
scalar component to the lighter component will be rapid, producing the lighter component and hard 
jets or leptons. The lighter component will again barrel through the detector producing a charged 
track and a late decay. This will also be a striking signal at the LHC. 

Since the x are always pair produced and decay at similar times, it may be possible to determine 
their lifetime by correlated measurements of double decay events in the calorimeter [61]. If all the 
particles in the multiplet are light enough to be produced at the LHC, a measurement of their mass 
and hfetime will be a direct probe of the GUT scale. The decay rate F of a particle of mass M 
through a dimension 5 operator scales as F ~ M^. A measurement of the masses and decay rates 
of the colored and electrically charged multiplets can then be used to determine if F and M scale 
as expected for a dimension 5 decay operator. Upon establishing this, a measurement of the scale 
A mediating the decay can be inferred from the relation F ~ 
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2. Fiveplet Relic 



The colored particles in (5, 5) will give rise to charged tracks and late time decays similar to the 
decays of the colored multiplets in (10, 10) as discussed earlier in subsection VII A 1. Electroweak 
symmetry breaking causes mass splittings between the electrically charged component /J and the 
neutral The collider phenomenology of the lepton doublet depends upon the spectrum of the 
theory. We first consider the case when the scalar components are lighter than their fermionic 
partners l^, so that the Z^, if produced, decay rapidly to the l^. Upon electroweak symmetry 
breaking, the masses of the scalar components Z+ and Z° are split by ~ cos (2/3) [55] . Depending 
upon the value of cos (2/3), the Z+ can be lighter than In this case, the Z+ will produce charged 
tracks as it traverses through the detector. Some of the Zjs will stop in the detector and result in 
late time decays. The will rapidly decay to the /° when the Z° is lighter than the l^. This decay 
produces hard jets and leptons along with missing energy from the /° that leave the detector and 
decay outside it. 

We now consider the case when the fermionic components are lighter than their scalar partners 
l^. In this case, the scalars l^, if produced, will rapidly decay to their fermionic partners l^. 
Electroweak symmetry breaking makes the electrically charged fermionic component heavier 
than the neutral component by ~ aMz ~ 350 MeV [53, 54]. The /^^ decays to the with a 
lifetime ~ (350^ McV)-''' ~ ^ x 10~^^ s = 1.2 cm yielding a displaced vertex. This displaced vertex 
produces soft pions, muons and electrons with energies ~ 350 MeV. The softness of these particles 
makes it difficult to trigger on them unless the production of these particles is accompanied by 
other hard signals like initial or final state radiation of photons [54]. One possible source for this 
hard signal is the decay of the scalars l-^. If the l-^ are produced, they will rapidly decay to their 
fermionic counterparts. The LSP will also be produced in this process resulting in the activation of 
missing energy triggers. 

Another possible source for this hard signal are the colored particles in the multiplet. In the 
model discussed in sub section VB, the colored particles in the 5*0(10) multiplet 16^^ can decay 
to the lepton doublets in the multiplet through dimension 5 operators. In this case, the late time 
decays of the colored multiplet can produce l^^ . This decay will be accompanied by a release of 
hadronic energy at the location of the colored particle and a charged track that follows the trajectory 
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of the produced This track ends in a displaced vertex when the decays to the The late 
time decay followed by the charged track can be used to trigger on this event. The discovery of 
soft particles at the location of the displaced vertex, while difficult experimentally, will unveil the 
presence of the lepton doublet. 

3. Singlet Relic 

The LHC signals for decays involving operators with singlet xs is dependent on the MSSM 
spectrum. If the MSSM LSP is heavier than x, then the MSSM LSP will decay to x at ~ 1000 
s. Cosmologically, this decay could solve the Lithium abundance problems and give rise to a dark 
matter abundance of x- Iii this scenario, since the MSSM LSP is no longer the dark matter, the 
MSSM LSP could be charged. The SUSY particles produced at the LHC will rapidly decay to the 
MSSM LSP and if this MSSM LSP is charged or colored, it will give rise to charged tracks and 
late decays in the detector. The signals of SUSY in this scenario are significantly altered since 
the decays of SUSY particles are no longer associated with missing energy signals as the charged 
MSSM LSP can be detected. However, if the MSSM LSP is neutral, it will decay to the x outside 
the detector and the LHC signals of this scenario will be identical to that of conventional MSSM 
models. 

The cosmological lithium problem could also be solved by decays of singlet xs to the MSSM if 
the X is heavier than the MSSM LSP. In this case, an abundance of x must be generated through non 
standard model processes. This abundance can be generated thermally through new interactions 
like a low energy U{1)b-l or through the decays of heavier, thermally produced standard model 
multiplets to x- In this scenario, these new particles could be discovered at the LHC, for example 
through Z' gauge boson searches. It is also possible that the initial abundance of x was generated 
through a tuning of the reheat temperature of the Universe. This scenario is devoid of new signals 
at the LHC. 
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X SU{h) Rep. 


Spectrum 


Prominent Signals and Features 


(10, 10) 


> Mlsp 


Charged tracks from colored and electrically charged particles. 
Measurement of mass, lifetime from stopped colored and charged particles. 
Potentially visible displaced vertex from colored doublet decays. 
Infer dimensionality and scale of decay from mass, lifetime measurements. 


(5,5) 


> Mlsp 


Charged tracks from colored particles. 

iVlCClo Lli ClliC-llLi Ui llidoo ^ 11 ic LllllC 11 Ulll o LUJJJJCVJ. t^L»lL»l CU. JJctl LICICO. 

Charged track from scalar lepton doublet. 

Potentially visible displaced vertex from lepton doublet decays. 


Singlet 


Mlsp > 


Possibility of charged or colored LSP giving rise to charged tracks. 
Charged LSP detection enables better measurement of SUSY spectrum. 




My, > Mlsp 


Discover new particles {e.g. TeV scale U{1)b-l) for thermally produced 



TABLE VIL Prominent signals of the dimension 5 decay operators considered in section IV. The signals 
are classified on the basis of the SU{5) representation of the new particle x its mass M^. When 
My. > Mlsp, the primordial Uthium problem is solved by the decays of x to the MSSM and vice- versa 
when Mlsp > My,. Bounds on dark matter direct detection [37] force % to decay to the MSSM LSP when 
X {e.g (5,5), (10, 10)) has dirac couplings to the Z. 



B. Signals of Dimension 6 Operators 

The dimension 6 GUT suppressed operators discussed in this paper mediate decays between 
singlet Ss and the MSSM when the operators conserve R-parity or cause decays of the LSP to 
the standard model when the operators violate R-parity (see section III) . The decay rate for these 
processes is F ~ 10~^^ s~^. The largest production channel for these particles at the LHC would 
be through the decays of colored particles. For example, the LSP would be produced in the decays 
of gluinos. The production cross-section for a 200 GeV gluino at the LHC is --^ 10^ fb[61]. With 
a integrated luminosity of 100 fb~^, the LHC will produce ~ 10^ LSPs through the decays of the 
gluino. This number is far too small to allow for the observation of the decay. Similarly, the number 
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of singlet Ss that can be produced either through a new low energy gauge symmetry {e.g. U (1) b-l) 
or through the decays of other particles that carry standard model quantum numbers is also too 
small to allow for the observation of the decay of the S. 

The most promising signatures of these dimension 6 GUT suppressed operators arc the astro- 
physical signatures of their decay. In section III, we show that these dimension 6 GUT suppressed 
decays can explain the observations of PAMELA/ ATIC. PAMELA has observed an excess in the 
positron channel while it does not see an excess in the anti proton channel. The hadronic branching 
fraction of these decays must be smaller than 10 percent in order to account for the PAMELA signal 
(see section III) . The need to suppress the hadronic branching fraction of these decays can be used 
to arrive at constraints on the supcrpartncr spectrum. 

For concreteness, wc consider the model in section III A where the scalar component s of the S 
gets a TeV scale vev (s). This leads to the operators and where s denotes the fermionic 

component of S, ^, ij represent slepton and lepton fields and (g, q) represent squark and quark 
fields. This operator will mediate the decay of s to the MSSM LSP if the mass Mg of s is greater 
than the LSP mass M^sp- The hadronic branching fraction of this decay depends upon the mass 
splittings AMgq and AM^^- between s and the squarks and sleptons respectively. The following 
options emerge: 

• Ms > Mf and Mg > Mf This decay produces on shell sleptons and squarks. The hadronic 
branching fraction is ~ ^ am" - ) ^^^^ fraction is smaller than 10 percent if AMgq < 
0.45 AM^i. 

• Ms < Ml and Mg < Mq-. This decay produces electrons and quarks through off-shell sleptons 
and squarks respectively. The hadronic branching fraction of this decay mode ij^j where 
Mf and Mq are slepton and squark masses respectively. This branching fraction is smaller 
than 10 percent if M|- < 1.8 Mq. 

• Ms > Mf and Mg < Mf This decay produces sleptons on shell. The dominant hadronic 
branching fraction of this decay mode arises from final state radiation of W and Z bosons 
yielding a branching fraction ~ 10^^ for Mg > 1 TeV. The hadronic branching fraction pro- 
duced as a result of decays through off-shell squarks is smaller than 10~^ since it is suppressed 
by additional phase space factors and couplings. 
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The cases where the slepton is hght enough to be produced on shell {e.g. Ms > Mf) allow for 
the possibility of another correlated measurement at the LHC and astrophysical experiments like 
PAMELA/ATIC. The decay of the s to a on-shell lepton and slepton will generate a primary source 
of hot electrons from the leptons produced in this primary decay. These electrons will produce 
a "bump" in the spectrum that will be cut off at roughly ~ The sleptons produced in this 
process are also unstable and will rapidly decay to a lepton and the LSP. This process produces a 
secondary lepton production channel that will also lead to a spectral feature cut off at M^— Mlsp- 
This astrophysical measurement can be correlated with measurements of this mass difference at the 
LHC and these independent measurements can serve as a cross-check for this scenario. 

The other possible decay in this model is the decay of the LSP to the s which will happen if 
Mlsp > Mg. In this case, leptons and quarks are produced through off-shell sleptons and squarks 
respectively (see figure 2) . The hadronic branching fraction of this decay is given by ij^j ■ This 
branching fraction is smaller than 10 percent if Mf < 1.8 Mq. We note that the model discussed 
in this section naturally allows for small hadronic branching fractions. The hadronic branching 
fraction is determined by slepton and squark masses and since squarks are generically expected 
to be heavier than sleptons due to RG running, the hadronic branching fraction of these models 
is generically small. These decays have a hadronic branching fraction that is at least ~ 10~^ for 
Ms ~ 1 TeV due to final state radiation of Z and W bosons. A measurement of anti-protons 
at PAMELA that indicates a hadronic branching fraction smaller than ~ 10~^ will rule out this 
model. A hadronic branching fraction that is larger than 10~^ will however constrain the squark 
and slepton masses as discussed earlier in this section. 



VIII. CONCLUSIONS 



In broad classes of theories, global symmetries that appear in nature, such as baryon number, 
are accidents of the low energy theory and are violated by GUT scale physics. If the symmetry 
stabilizing a particle, such as the proton, is broken by short-distance physics then it will decay with 
a long lifetime. In a sense, a particle is continually probing the physics that allows it to decay. 
Thus, though the effects of GUT scale physics are suppressed at low energies, this suppression 
is compensated for by the long time scales involved. Intriguingly, short-distance physics can be 
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unveiled by experiments sensitive to long lifetimes. 

Astrophysics and cosmology provide natural probes of long lifetimes. The dark matter particle 
appears stable, but dimension 6 GUT suppressed operators can cause it to decay with a long lifetime 
~ 10^^ s. This can lead to observable signals in a variety of current experiments including PAMELA, 
ATIC, HESS, and Fermi. Such observations can also probe the TeV scale physics associated with 
dark matter annihilations. However, annihilations cannot proceed through GUT scale particles 
since such a cross section would be far too small. Thus, the observation of annihilations can probe 
TeV physics but does not probe GUT physics directly. 

This leads to interesting differences between the expected signals from dark matter decays 
and annihilations [53, 72]. Decaying dark matter can directly produce photons with an 0{1) 
branching ratio, as seen in Section HI. It can also produce light leptons without helicity or p-wave 
suppression and give signals in the range necessary to explain the PAMELA/ ATIC excess without 
boost factors [35, 73, 74, 75, 76]. Additionally, in our framework the hadronic branching fractions 
can be naturally suppressed since sleptons tend to be lighter than squarks, perhaps explaining 
the lack of an antiproton signal at PAMELA. We have also found that the signals produced by 
decays may show many interesting features beyond just a bump in the spectrum of electrons and 
positrons. For example, the spectrum from decays may accommodate the secondary feature visible 
in the ATIC data at energies between 100 - 300 GeV. The SUSY spectrum can allow for the direct 
production of a slepton, lepton pair from the decaying dark matter. In this case, the subsequent 
decay of the slepton provides a secondary hot lepton that can give rise to the observed secondary 
ATIC feature. It is also common in our framework to have two dark matter particles, scalar and 
fermion supcrpartners, that both decay at late times producing interesting features in the observed 
spectra of electrons or photons. 

Such models will be tested at several upcoming experiments. The spectrum of electrons and 
positrons will be measured to increasing accuracy by Fermi, HESS, and PAMELA, testing the 
observed excesses. These measurements could potentially reveal not just the nature of dark matter, 
but also the mass spectrum of supersymmetry and GUT scale physics. Excitingly, Fermi and 
ground-based telescopes such as HESS or MAGIC will be providing measurements of the gamma- 
ray spectrum in the near future. These measurements have the potential to test many of our 
models, and in particular may be able to observe final state radiation from models that explain the 
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electron/positron excess. Further, they could distinguish decaying and annihilating dark matter 
scenarios by the different angular dependence of the gamma ray signals. 

The supersymmetric GUT theories in which the dark matter decays through dimension 6 oper- 
ators often have TeV mass particles that decay via dimension 5 operators with a lifetime of ~ 1000 
sec. This results in relic decays during BBN which could explain the observed primordial Lithium 
abundances. This TeV mass particle can be produced at the LHC and if charged or colored will stop 
and decay inside the detector within ~ 1000 sec. This will allow the measurement of its properties 
and directly make the connection with primordial nucleosynthesis. 
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